CHARACTERIZATION OF THE NEWTONIAN 
FREE PARTICLE SYSTEM 
IN m ^ 2 DEPENDENT VARIABLES 



JOEL MERKER 

Abstract. We treat the problem of linearizability of a system of second order or- 
dinary differential equations. The criterion we provide has applications to nonlinear 
Newtonian mechanics, especially in three-dimensional space. 

Let K = M or C, let x G K, let m ^ 2, let y := (y^, . . . , y") G K'" and let 

vlx = {x,y,yx) , , vTx = F"' i^, y, Vx) , 

be a collection of m analytic second order ordinary differential equations, in general 
nonlinear. We obtain a new and applicable necessary and sufficient condition in order 
that this system is equivalent, under a point transformation 

(x, y\ . . . , y^) ^ {X{x, y), Y\x,y),..., y™(x, y)) , 

to the Newtonian free particle system Y^j^ — ■ ■ ■ — YJ^-^ — 0. 

Strikingly, the explicit differential system that we obtain is of first order in the case 
TO ^ 2, whereas according to a classical result due to Lie, it is of second order the 
case of a single equation (m = 1). 
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§1. Introduction 

Several physically meaningful systems of ordinary differential equations are of sec- 
ond order, as for instance the free particle in m-dimensional space, the damped or 
undamped harmonic oscillator, coupled or not, having constant time-dependent fre- 
quency or not, etc. Such systems are ubiquitous in Newtonian Mechanics, in Hamil- 
tonian Dynamics and in General Relativity. 

Two classical major problems are to classify these systems up to point or contact 
equivalence (Lie's Grail) and to recognize when they coincide with the Euler equa- 
tions associated to a Lagrangian (inverse variational problem). In small dimensions, 
complete results hold (Lie, Tresse, Cartan; Darboux, Douglas). However, in arbitrary 
dimension, both tasks quickly exceed the human as well as the digital computer scale. 
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due to the intrinsic complexity of the underlying symbolic computations (explosion, 
swelling) and to the exponentially increasing number of cases to be treated. We refer 
to Olver's monograph IIOL1995I I for a panorama of problems, methods and results. 

At least, as a first step in classification, with respect to applications, there are both 
a mathematical and a physical interest in determining concrete, explicit and applicable 
("ready-made") criteria for a system of ordinary differential equation to be equivalent, 
via a local point transformation, to a linear equation. 

1.1. Scalar equation. In this respect, we remind the celebrated linearizability criterion 
for a single equation, due to Lie. Let K = M of C. Let a; G K and y G K. Consider a 
local second order ordinary differential equation y^x = F{x, y, yx), possibly nonlinear, 
with a locally K-analytic right-hand sidqj. 

Theorem 1.2. ( [|Liel883L pp. 362-365; IIGTW1989II : I10L1995L p. 406) The follow- 
ing four conditions are equivalent: 

(1) y-xx = F{x,y,yx) is equivalent under a local point transformation {x,y) ^ 
{X, Y) to the free particle equation Yxx = 0; 

(2) yxx = F{x,y,yx) is equivalent to some linear equation Yxx = Go{X) + 
Gi{X)Y + H{X)Yx; 

(3) the local Lie symmetry group of yxx = F{x, y, yx) is eight-dimensional, locally 
isomorphic to the group PGL(3, K) of all projective transformations ofP2(K). 

(4) Fy^ya^y^ya^ = 0, or equivalently F = G + yxH + {yxf L + {yxf M, where 
G, H, L, M are functions of{x, y) that satisfy: 

^ = -2Gyy + ^ Hxy -^Lxx + 2{GL)y-2GxM -AGMx + ^HLx- ^HHy, 
^ = -\Hyy + \Lxy-2Mxx + 2GMy+AGyM -2{HM)x-\HyL+\LLx. 

Section 2 of this paper is devoted to a detailed exposition of the original proof 
of the equivalence between (1) and (4), following [|Liel883l . In the contemporary 
literature, to the author's knowledge, there is no modem restitution of Lie's elegant 
proof, whereas the description of an alternative proof of Theorem 1.2 as a byproduct 
of E. Cartan's equivalence algorithm appears in the references nTR1896ll . HCa 192411 , 
[IGTW1989II . [IOL1995L IINS2003fl 

Lie's Grail would comprise: 

• a complete classification of all Lie algebras of local vector fields; Lie achieved 
this task in dimension 2 over C (real case: HGKO 199211 ); however, as soon 

^By borrowing techniques developed in IIMa2003l . this theorem as well as the next both hold under 
weaker smoothness assumptions, namely with a or a W^^^ right-hand side. 

We note that in these references, the already substantial computations are stopped just after the re- 
duction to an {e}-structure on an eight-dimensional (local) principal bundle over the three-dimensional 
first order jet space. The vanishing of two (among four) fundamental tensors in the structure equations 
of the obtained {e}-structure yields two partial differential equations satisfied by the right-hand side 
F{x, y, Ux), which are equivalent to (4) of Theorem 1.2. We mention that with the help of Maple pro- 
gramming, the complete reduction to an {e}-structure on the base (not only on the principal bundle) is 
achieved in IIHK1989L in the simpler case of so-called fiber-preserving transformations, namely point 
transformations leaving invariant the "vertical" foliation {x = ct.}. To the author's knowledge, the 
complete confirmation of Tresse's results by means of E. Cartan's method has never been achieved. 
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as the dimension is ^ 3, the complete classification is unknown, due to the 
intrinsic richness of imprimitive Lie algebras of vector fields; 

• a list of all the possible Lie algebras that can be realized as infinitesimal Lie 
symmetry algebras of partial differential equations, together with their Levi- 
Malcev decomposition; 

• an explicit Grobner basis of the (noncommutative) algebra of all differential 
invariants of each equation in the list. 

However, complete results hold only for the scalar second order ordinary differen- 
tial equation. Some tables extracted from Lie's Gesammelte Abhandlungen and from 
Tresse's prized thesis UTR 189611 may be found in IIOL1995L 

1.3. Systems. Let x G K, let m ^ 2 , let y := (y^ . . . , y™) e K*" and let 

(1.4) yUx) = F' (x, y{x), y,{x)) , , yZ{x) = (x, y{x), y,{x)) 

be a collection of m analytic second order ordinary differential equations, possibly 
nonlinear, of the most general form. 

Motivated by physical applications and by geometrical questions, some au- 
thors (Leach [|Lel980L Grissom-Thompson-Wilkens IIGTW1989L Gonzalez- 
Lopez I1GL1988II Pels [|Fel995ll, Crampin-Martm ez-Sarlet I1CMS1996L 
Doubrov [|Do2000L Grossman [IGr2000L Mahomed-Soh IIMS2001L and others) 
have been interested in a least characterizing those that have the Lie symmetry group 
of maximal dimension. In llLel980l , based on the belief that the equivalence between 
(1), (2) and (3) of Theorem 1.2 would persist in the case m ^ 2, it was conjectured 
that the symmetry algebra of every linear system 

m m 

(1.5) yi, = Gi{x) + y'G{^,{x) + y'^ H^^{x), j = l,...,m, 

1=1 h=i 

of m ^ 2 second order differential equations has a Lie symmetry group locally isomor- 
phic to the full transformation group PGL(m + 1, K) of the projective space Pm+i(K). 
However, Gonzalez-Lopez I1GL1988I infirmed this expectation, and produced a neces- 
sary and sufficient condition (see Corollary 1.8 below) for local equivalence of such 
linear systems to the free particle system 

(1.6) >lx = 0, j = l,...,m. 

Applying Lie's algorithm it is easily seen ( [IGG1983I ) that the free particle system has 
a local symmetry algebra of dimension < + 4 m + 3, the bound being attained by 
y^^ = 0, J = 1, . . . , m, with group PGL(m + 1, K). 

In 1939, for fiber-preserving transformations only, Chern llChl939ll conducted the 
E. Cartan algorithm through absorptions of torsion, normalizations and prolongations 
up to the reduction to an {e}-structure. Remarkably, in 1995, Pels llFel995ll conducted 
the E. Cartan algorithm for general systems y^^ = F^{x, y, yx) and for general point 
transformations. As a byproduct of the uniqueness of the obtained {e}-structure for 
which all invariant tensors vanish. Eels deduced in nEel995ll that the flat system Y^^^ = 
0, j = 1, . . . , m, is, up to equivalence, the only system of second order possessing a 
symmetry group of maximal dimension. Alas, the {e}-structures obtained by Chem 
and by Eels are not parametric, so that the counterpart to (4) of Theorem 1.2 was 
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lacking as soon as m ^ 2. The extraordinary heaviness of E. Cartan's method is a well 
known obstacle. 

Recently, Neut IIN2003II (see also [IDNP2005II ) wrote a Maple program to compute 
(among other {e}-structures) Pels' tensors in a parametric, explicit way. The digital 
computations succeeded in the case m = 2, yielding Theorem 3 of [IDNP2005II . a 
statement that may be checked to be equivalent to Theorem 1.7 (3) just below in the 
case m = 2. In the physically most meaningful case m = 3, a present-day computer is 
stuck ( I1N2003II . I1DNP2005II ). 

Extending Lie's less heavy computations El, we present here a complete solution to 
the characterization of Y^^ = for arbitrary m ^ 2. 

Theorem 1.7. Suppose m ^ 2. The following three conditions are equivalent: 

(1) the system yl^ = F^{x, y, yx), j = 1, . . . ,m, is equivalent, under a local point 
transformation (x, y^) t— > (X, Y^ ) to the free particle equation Y^-^ = 0; 

(2) the local Lie symmetry group of yl^ = {x,y,yx) is {rn? + 4 m + 3)- 
dimensional and locally isomorphic to PGL(m + 1, IK); 

(3) the right hand sides F^{x, y, yx) are of a special form, described as follows. 

(i) There exist local "K- analytic functions G\ Hj^, Lj_^ and Mi^^i^, where 
j,h,h = I, ... ,171, enjoying the symmetries Lj^ = Lj^ and Mi-^^i^ = 
Mi^^i^ and depending only on (x, y) such that F^{x, y, yx) may be written 
as the following specific cubic polynomial with respect to yx'. 

m mm mm 

/l=l h = l l2 = l h = l l2 = l 

(ii) The functions , H^^, <^nd Mi^^i^ satisfy the following explicit system 
o/four families of first order partial differential equations: 

{ = - 2G;,, +251 G% + Hl^, - 51 H\1,+ 



(I) 



m 



+ 2j2G'Ll,-2 5lJ2G'Ltk+ 

k=l k=l 
^ m ^ m 

k=l k=l 



Throughout the article, we do not adopt the summation convention, because in several subsequent 
equations, some repeated indices shall appear that will not be summed. Also, we always put commas 
between the indices. For instance ^ denotes dLj^ i^/dy^^ shortly. As usual, 51 is the Kronecker 
symbol. 



(11) 



(III) 



(IV) 
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where the indices j, h vary in {1,2, ... , m}; 







-^ — 



1 



1 



I . ^ 2-1 



51 Ml 

k=l k=l 

^ m ^ m 

k=l k=l 

\ k=l k=l / 

/ m m \ 



where the indices j, h, I2 vary in {1,2, ... , m}; 



k=l k=l 

^ m ^ m 

+ -2<Y. Hi - 2 E Hi M,3,,+ 

k=l k=l 

m m 

fe=l fc=l 

where the indices j, li, I2, 13 vary in {1, . . . m}; and 



= M, 



M, 



h,h,y''2 



k=l 



k=l 



where the indices h,l2, h vary in {1, ... , m}. 



Let us provide commentaries and explanations. The form of the right-hand side of 
(3)(i) is the analog of the form of the right-hand side F in (4) of Lie's Theorem 1.2. 
However, we notice that the right-hand side of (3)(i) is not the most general degree 
three polynomial in the variables yl, j = 1, . . . ,m: some coefficients of the cubic 
terms vanish. 

Very strikingly, the differential system (I), (II), (III), (IV) satisfied by the functions 
G^, Hl^, Mi^^i^ is oi first order for m ^ 2, whereas the system (4) satisfied by 

G, H, L, M in Lie's Theorem 1.2 is of second order for m = 1. This confirms the 
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main theorem of Gonzalez-Lopez IIGL1988I I. that is recovered here in the special linear 
case (1.5). 

Corollary 1.8. ( [IGL1988I ) For m ^ 2, a linear (nonhomogeneous) system yl^ = 
Gq{x) + G{ i{x) + J2h=i -^/i(^) equivalent to Y^-^ = if and only if 

there exists a function B{x) such that the m x m matrix Gi may be written under the 
specific form: 

1 1 " 

(1-9) G\, = -Hi--Y.HiHt + 5lB. 

k=l 

We also mention that Fels obtained a characterization of equivalence to the free par- 
ticle Fj^jf = 0, j = 1, . . . , m, by the vanishing of two (nonexplicit) tensors Sfj^i and 
P- (Corollary 5.1 in [|Fel995ll ). In this reference, some parametric computations are 
achieved after restricting the initial G-structure, together with its subsequent prolonga- 
tions, to the identity element of the group; explicit expressions of Sjf^i]-^^ and of -P/|jj 
are then obtained, through already hard computations. The vanishing of the two tensors 
P- and S^j^i at the identity of the structure group, namely, as computed in Lemma 4.1 
of l|Fel995l , yields (translating into our notation) 
(1.10) 



Id VxVxyx n + 2 ^-^ ^-^ ^ ' vxVx vx 



ll a-(i) CT(fc)5 



0= {Pi 



Id 




4 2^ 2^ ^yi ^y'l 



k=l 1=1 

where D is the total differentiation operator Ux ^~^Ym^i ^^'^ where 

i, j, k, I = 1, . . . , m. Strikingly, one may check that the first equation is equivalent 
to (3)(i) of Theorem 1.2 and then that the second equation yields the (complicated) 
four families of first order partial differential equations (I), (II), (III) and (IV). So, in 
Corollary 5.1 of l|Fel995l , one may replace the vanishing of S^-j^^ and of P- by the 
vanishing of S^-i^i\^^ and of -P/ [j^, which were explicitly computed there. 

This phenomenon could be explained as follows: as soon as the tensors Sjj^i vanish, 
the system enjoys a projective connection (appendix of [|Fel995l ): with such a connec- 
tion, the tensors P- then transform according to a specific rule via tensorial rotation 
formulas and their general expression may be deduced from their expression at the 
identitjQ. We have checked this, but as we try to avoid the method of equivalence, de- 
tails will not be reproduced here. Similar observations appear in Hachtroudi [|Hal937l . 

Even if the expressions (1.10) are more compact than the (equivalent) conditions in 
Theorem 1 .7 (3), we prefer the complete expressions of Theorem 1 .7 (3), since they 
are more explicit and ready-made for checking whether a physically given nonlinear 



^ Similar rotation formulas are known in the much simpler case of (pseudo-) Riemannian metrics, 
see Chapter 12 of Olver MOL 19951 . It would be interesting to write a program, finer and more effi- 
cient than IIN2003I . which would systematically recognize such rotation formulas in any application of 
E. Cartan's equivalence method. 
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system is equivalent to a free particle. If the reader prefers compact expressions and 
"short" theorems, (s)he may replace the conditions of Theorem 1.7 (3) by (1.10). 

Open Problem 1.11. Characterize explicitly the linearizability of a Newtonian system 
in m ^ 2 degrees of freedom, i.e. local equivalence to : 

m m. 

(1.12) = Gi{X) + Y^G{,{X) + Yx Hl{X)- 

1=1 h=i 

1.13. Organization, avertissement and acknowledgment. Section 2 is devoted to 
a thorough restitution of Lie's original proof of the equivalence between (1) and (4) 
in Theorem 1.2. Section 3 is devoted to the formulation of combinatorial formulas 
yielding the general form of a system equivalent to Y^-^ = 0, j = 1, . . . , m, under a 
local K-analytic point transformation (x, y^) t— > {X, Y^), for general m ^ 2 ; the proof 
of the main technical Lemma 3.32 is exposed in Section 5. Section 4 is devoted to the 
final proof of the equivalence between (1) and (3) in Theorem 1.7, the equivalence 
between (1) and (2) being already proved by Fels [|Fel995ll . 

Some word about style and intentions. We wanted the proof of the equivalence 
between (1) and (3) be totally complete, every tiny detail being rigorously and patiently 
checked. This is why we decided to carefully detail each intermediate computational 
step, seeking first the combinatorics of the formal calculations in the case m = 1 and 
devising then the underlying combinatorics for the case m ^ 2. Actually, the size of 
differential expressions is relatively impressive, as will become soon evident. Thus, 
no intermediate symbolic computation will be hidden, hence essentially no checking 
work is left to the reader, as would have been the case if we did not have typed all the 
computations. 

We also would like to point out that except in every specific standard situations, as 
with the much studied Riemann and Ricci tensors, present-day computer programs are 
not yet powerful enough to apply the E. Cartan equivalence method when the num- 
ber of some collection of variables is a general integer. All the formulas obtained in 
Sections 2, 3 and 4 were first treated completely by hand and then, some of them were 
confirmed afterwards with the help of MAPLE release 6 in the cases m = 2 and m = 3. 
The author is indebted to Sylvain Neut and to Michel Petitot, from the University of 
Lille 1, for their help in computer machine confirmations. 

§2. Proof of Lie's theorem 

2.1. Argument. This preliminary section contains a detailed exposition of Lie's orig- 
inal proof of the equivalence between (1) and (4) in Theorem 1.2. Since our goal is to 
guess the combinatorics of computations in several variables, it will be a crucial point 
for us to explain thoroughly and patiently each step of Lie's computation. Without 
such an intuitive control, it would be hopeless to conduct any generalization to several 
variables. Hence we shall respect a fundamental principle: always explain clearly and 
completely what sort of computation is achieved at each step. Also, we shall many 
times introduce some appropriate new notation. 

2.2. Combinatorics of the second order prolongation of a point transformation. 

Let K = M or C. Let (x, y) i— > {X(x,y),Y(x,y)) be a local K-analytic invertible 
transformation, defined in a neighborhood of the origin in K^, which maps the sec- 
ond order differential equation y^x = F{x,y,yx) to the fiat equation Yxx = 0. By 
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assumption, the Jacobian determinant 
(2.3) A{x\y) :-- 



Xx Xy 

Y Y 



is nowhere vanishing. Since the equation Yxx = is left unchanged by any affine 
transformation in the (X, Y) space, we can (and we shall) assume that the transforma- 
tion is tangent to the identity at the origin, namely the above Jacobian matrix equals 
the identity matrix at {x,y) = (0,0). 

The computation how the differential equation in the (X, Y) coordinates is related to 
the differential equation in the (x, y) -coordinates is classical, cf. llLiel883L I1TR1896L 
I1BK1989L HIB 199211 : let us remind it. A local graph {y = y{x)} being transformed to 
a local graph {Y = Y{X)}, we have a direct formula for the first derivative Yx'. 

dY dx- dY{x, y{x))/dx Y^ + yxYy 

(2.4) Yx ~ ~ 



dX dx ■ dX{x,y{x))/dx X^ + yxXy 

This yields the prolongation of the transformation to the first order jet space. For the 
second order prolongation, introducing the second order total differentiation operator 
(which geometrically corresponds to differentiation along graphs {{x, y{x))}) defined 
by 

^^ ^ d d d 

(2.5) D:= — + yx — + y. 



OX oy oyx 

we may compute, simplify and reorder the expression of the second order derivative in 

the (X, F) -coordinates: 

(2.6) 

' d^Y _DYx _D [{Yx + yxYy){Xx + yxXx)-'] _ 



Y, 



XX ■- 



dX^ DX Xx + yxXy 

7T7 I Y' 13 {y^^ \XxYy YxXy^ -\- XxYxx YxXxx~\~ 
[■^x I Vx^yl 

+ Ux ['^{XxYxy — YxXxy) — {XxxYy — YxxXy)] + 
VxVx \^xYyy YxXyy '2{XxyYy Y^j^Xj^)] -|- 

+ yxyXyx[-{XyyYy -YyyXy)]} . 

Even if not too complicated, the internal combinatorics of this expression has to be 
analyzed and expressed thoroughly. First of all, as Yxx = by assumption, we may 
erase the cubic factor [Xx + yxXy]~^. Next, as the factor of yxx in the right-hand side 
of (2.6), we just recognize the Jacobian A(x|i/) expressed in (2.3) above. Also, all the 
other factors are modifications of the Jacobian A(x|y), whose combinatorics may be 
understood as follows. 

There exist exactly three possible distinct second order derivatives: xx, xy and yy. 
There are also exactly two columns in (2.3). By replacing each of the two columns of 
first order derivative in A(a;|y) by any column of second order derivative (leaving X 
and Y unchanged), we may build exactly six new determinants 



(2.7) 



/\{xx[y) /^{xy[y) A{yy[y) 
lS.{x[xx) A(x|xy) A{x[yy) 
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where for instance 
(2.8) ( ^{xx\y) :-- 



Y Y 



and A(x|x?/) : = 



X, 
Y 



xy 



xy 



Hence, by rewriting (2.6), we see that the equation y^x = F{x, y, y^) equivalent to 
Yxx = may be written under the general explicit form, involving determinants 
(2.9) 







yx 



Xx 



+ Vxyx 



Xy 

X, 



Xx 



Xx 

Y 



X, 



yy 



yy 



Vx 



Y Y 

-^xy y 

Y Y 

J-xy J-y 



Xx 



X, 
Y. 



xy 



xy 



+ VxVxVx 





Y 

-^xx 


Xy 


1 




Y 

^ XX 


Yy 




L 


Y 

^yy 


Xy 


1 




Y 

^yy 


Yy 





or equivalently, after solving in yxx, i-S- after dividing by the Jacobian A(a;|?/): 



(2.10) 



A(x|xx) 
~ A(x|y) 

+ {yx)' 



yx 

A(x|yy) 



A(x|xy) A(xa;| 



A(x|y) A(x|y) 



+ 



^{yy\y) \ 

A(x|y) / 



At this point, it will be convenient to slightly contract the notation by introducing a 
new family of square functions as follows. We first index the coordinates {x,y) as 
{y^,y^), namely we introduce the two notational equivalences 



(2.11) 



y = X, y' 



y 



which will be very convenient in the sequel, especially to write down general combi- 
natorial formulas anticipating our treatment of the case of m ^ 2 dependent variables 



hand, our six square functions U'^]^ , symmetric with respect to the lower indices, 



to be achieved in Sections 3, 4 and 5 below. With this convention at 

-ifci 



where ^ ji,j2,ki ^ 1, are defined by 

A{xx\y) 



(2.12) 



□0 := 

— XX 



A{x\y) 
A{x\xx) 
A(x\y) 



□0 

'—'xy 



xy 



M^y\y) 

A{x\y) 
A{x\xy) 
A(x|y) ' 



a 



yy ■ 



yy 



A{x\y) ^ 
A(x|yy) 
A{x\y) ■ 



Here of course, the upper index designates the column upon which the second or- 
der derivative appears, itself being encoded by the two lower indices. Even if this 
is hidden in the notation, we shall remember that the square functions are explicit 
rational expressions in terms of the second order jet of the transformation (x, y) ^ 
{X{x,y),Y{x,y)). However, we shall be aware of not confusing the index in the 
square functions with a second order partial derivative of some function denoted 
by the square symbol: indeed, the partial derivatives are hidden in some determinant. 
At this point, we may summarize what we have established so far. 

Lemma 2.13. The equation yxx = F{x,y,yx) is equivalent to the flat equation Yxx = 
if and only if there exist two local K- analytic functions X{x, y) and y) such that 
it may be written under the form 

(2.14) y^^ = -Ul^ + yx-{- 



-2 nly + DL) + {y^f ■ {-^ly + 2 J + {yxf ■ □: 



yy 
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At this point, for heuristic reasons, it may be useful to compare the right-hand side 
of (2.14) with the classical expression of the prolongation to the second order jet space 
of a general vector field of the form L := X{x,y)-^ + Y(x,y)-^, which is given, 
according to llLiel883L II011986II . I1BK1989L by ^ 
(2.15) 

L^'^ =X^ + Y^+[Y^ + y^.iYy-X^) + {y^f ■ (-X,)] A+ 

d 

+ Vxx ■ {Yy - 2X^) + y^y^^ ■ {-3Xy)] . 

(yyxx 

We immediately see that (up to an overall minus sign) the right-hand side of (2.14) 
is formally analogous to the second line of (2.15) : the letter X corresponds to the 
symbol and the letter Y corresponds to the symbol This analogy is no mystery, 
just because the formula for L*^^^ is classically obtained by differentiating at e = the 
second order prolongation [exp(£:L) (■)]^^^ of the flow of L ! 

In fact, as we assumed that the transformation {x, y) t— > {X{x^ y), Y{x, y)) is tan- 
gent to the identity at the origin, we may think that X^ = 1, Xy = 0, Y^ = and 
Yy = 1, whence the Jacobian A{x\y) = 1 and moreover 



(2.16) 




|— |U Y 
yy ~ ^yy^ 

— Y 
^yy ~ ^yy 



By means of this (abusive) notational correspondence, we see that, up to an overall mi- 
nus sign, the right-hand side of (2.14) transforms precisely to the second line of (2.15). 
This analogy will be useful in devising combinatorial formulas for the generalization 
of Lemma 2. 13 to the case of m ^ 2 variables (y^, . . . ,y"^), see Lemmas 3.22 and 3.32 
below. 



2.17. Continuation. Clearly, since the right-hand side of (2.14) is a polynomial of 
degree three in y^, the first condition of Theorem 1.2 (4) immediately holds. We are 
therefore led to establish that the second condition is necessary and sufficient in order 
that there exist two local K-analytic functions X(x,y) and Y{x,y) which solve the 
following system of nonlinear second order partial differential equations (remind that 
the second order jet of (X, Y) is hidden in the square functions): 



( G = - Dl 



(2.18) 



H= - 



2 + 



^^xy^ 



M 



yy 



In the remainder of this section, following [|Liel883l l. p. 364, we shall study this sec- 
ond order system by introducing two auxiliary systems of partial differential equations 
which are complete, and we shall see in §2.38 below that the compatibility conditions 
(insuring involutivity, hence complete integrability) of the second auxiliary system ex- 
actly provide the two partial differential equations appearing in Theorem 1.2 (4). 
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2.19. First auxiliary system. We notice that in (2.18), there are two more square 
functions D^.^, D^y, Hyy, H].^, nl.y, Hyy, than functions G, H, L and M. Hence, 
as a trick, let us introduce six new independent functions n^,^ ,,2 (•^' symmetric 
with respect to the lower indices, for ^ ji, j2, ^1 ^ 1 and let us seek necessary 
and sufficient conditions in order that there exist solutions {X, Y) to the first auxiliary 
system: 

-'-'■0,0' '-'xy ~ ■'■■'■0,1' '-'yy ~ ■'■■'■1,1' 

ttI — TT^ — TT^ 

■'■■'■0,0' '-'xy — ■'■■'■0,1' '-'yy ~ ■'■■'■1,1- 

According to the (aprooximate) identities (2.16), this system looks like a complete 
second order system of partial differential equations in two variables (x, y) and in 
two unknowns (X, Y). More rigorously, by means of elementary algebraic operations, 
taking account of the fact that X^ = 1, Xj, = 0, = and Yy = 1, one may transform 
this sytem in a true second order complete system, solved with respect to the top order 
derivatives, namely of the form 

ao V — \^ V — \^ 

^^■0,0' ^xy — ^io,i' ^yy ~ ^^■1,1) 

^^■0,0' ^xy ~ ^^^0,1' ^yy ~ ^^■i,i) 

where the A^^^ are local K-analytic functions of (x, y, X, Y, X^, Xy, Y^, Yy). For such 
a system, the compatibility conditions [which are necessary and sufficient for the exis- 
tence of a solution {X, Y)] are easily formulated: 

(AJ,,)., (AJJ, = {Al 

Equivalently, we may express the compatibility conditions directly with the sys- 
tem (2.20), without transforming it to the form (2.21). This direct strategy will be 
more appropriate. 



(2.20) 



— XX 



(2.21) 



Xt 



Y 



(2.22) 



(^o,o)y 
(^o,o)y 



2.23. Compatibility conditions for the first auxiliary system. Indeed, to begin with, 
let us remind that the A(-|-) are determinant, hence we have the skew-symmetry re- 
lation A(x"'y''\x'^y^) = —A(x'^y'^\x"'y'') and the following two formulas for partial 
differentiation 

r rA(xVkV)l = A(x'^+^/|xV) + A(xV|x^+^/), 

(2 24) < 

[ [A(xV|xV)]^ = A(xV+W) + A(xV|xV+^). 

With these formal rules at hand, as an exercise, let us compute for instance the follow- 
ing cross differentiation (remember that the lower index in the square functions is not 
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a partial derivative): 



(2.25) 



CnM - fnM = — ( ^^^^^ ^My) \ ^ 



[A(x 



[A(x|y)]^ 



A(xxy|y) ■ A(x|y) g+ /^{xx\yy) ■ A(x|?/)- 

- A(xy|y) ■ A(xx||/) ^- A(x|y?/) ■ A(xx|?/)- 
- A(3:xy|y) ■ A(x|?/) q- A(a:|/|x|/) ■ A(x||/) q+ 

+ /\{xy\y) ■ A(xx||/) ^+ A(x?/|y) ■ A(x|x?/) 

{A(xx|?/2/) ■ A(a;|2/) - ^{x\yy) ■ ^{xx\y) + 
+/^{xy\y) ■ A(x|x?/)} . 



Crucially, we observe that the third order derivatives kill each other and disappear, see 
the underlined terms with ©appended. Also, two products of two determinants A(-|-) 
involving a second order derivative upon one column of each determinant kill each 
other: they are underlined with (gi appended. Finally, by antisymmetry of determinants, 
the term ^{xy\xy) ■ A(a;|y) vanishes gratuitously: it is underlined with ©appended. 

However, there still remains one term involving second order derivatives upon the 
two columns of a determinant: it is A(xx\yy). 

We must transform this unpleasant term A{xx\yy) ■ A(x\y) and express it as a prod- 
uct of two determinants, each involving a second order derivative only in one column. 
To this aim, we have: 

Lemma 2.26. The following three relations between the differential determinants 
A(-|-) hold true: 



{121) 



A.{xx\xy) ■ A(x|?/) = /S.{xx\y) • A(x|xy) 
A{xx\yy) ■ A{x\y) = A{xx\y) ■ A{x\yy) 
A{xy\yy) ■ A{x\y) = A{xy\y) ■ A{x\yy) 



A{xy\y) ■ A{x\xx), 

Hyy\y) ■ ^ix\xx), 
^{yy\y) ■ ^{x\xy). 



Proof. Each of these three formal identities is an immediate direct consequence of the 
following Pliicker type identity, easily verified by developing all the determinants : 
(2.28) 





Bi 




Ci 


Di 




Ai 


Di 




Ci 


Bi 




Bi 


Di 




Ci 


Ai 


A2 


B2 




C2 


D2 




A2 


D2 




C2 


B2 




B2 


D2 




C2 


A2 



where the variables Ai, A2, Bi, B2, Ci, C2, -Di, -D2 £ IK are arbitrary. 



□ 



Thanks to the second identity (2.27), we may therefore transform the result left 
above in the last two lines of (2.25); as desired, it will remain determinants having 
only one second order derivative per column, so that after division by [A(a;|y)]^, we 
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discover a quadratic expression involving only the square functions themselves: 
(2.29) 

(□1)^ - = {^{xx\y) ■ ^{x\yy) - ^{yy\y) ■ ^{x\xx)- 

-A(x|yy) ■ A(xx|?/) + A(xy|?/) ■ /^{x\xy)] 

= fa/i m2 {-^{yy\y) ■ + /^{xy\y) ■ A(x|x?/)} 

[A(x 



= -0° ■ + 0° ■ . 

— ^?/J/ — ^2;a; ' — xy — xy 

In sum, the result of the cross differentiation {^^xx)y ^ is a quadratic expression 

in terms of the square functions themselves! Following the same recipe (with no sur- 
prise), one may establish the following relations, listing all the compatibility conditions 
(the first one is nothing else than (2.29)): 



(2.30) 



Instead of checking patiently each of the remaining three cross above cross differenti- 
ation identities, it is better to establish directly the following general relation. 

Lemma 2.31. Remind from (2.11) that we identify y^ with x and y^ with y and let 
^ ji,j2,j3,fci ^ 1- Then 

1 1 
(2.32) (u\ A =-y U\ ^.-rt + V U\ ^-U^ 



r my - 


Hy)x = 


— XX 


yy 




^xy 










Hy)y- 


Hy)x = 


^xy 


yx 


-^ly 


°ly 


^ yy 


□L 




^ly 




^x = 


— XX 




-□L 


■°ly 


+ n° 


□L 


+ °ly 


■^ly 


.Hy)y- 


Hy)x = 


3;y 




+°;. 


— XX 











" k2=0 k2=0 



uhyH ■ 



This lemma is left to the reader; anyway, we shall complete the proof of a gener- 
alization of Lemma 2.31 to the case of m ^ 1 dependent variables (y^, . . . ,?/'") in 
Section 2 below (Lemma 3.40). 

Coming back to the first auxiliary system (2.20), we therefore have obtained a nec- 
essary and sufficient condition for the existence of {X, Y): the functions II^i p should 
satisfy the following system of first order partial differential equations, just obtained 
from (2.30) by replacing the square functions by the Pi functions: 
(2.33) 

-■-^o.o -'-'■1,1 ^ -'-'0,1 -'-'-0,1' 

^-'-0,1 -'-'0,1 -'-'-0,1 -'-'-1,1 ^ -'-'-1,1 -'-'-0,0 ^ -'-'1,1 -'-'-0,1! 
-'-'-0,0 -'-'0,1 -'-'-0,0 -'-'-1,1 ^ -'-'-0,1 -'-'1,1 ^ -'-'-0,1 -'-'-0,1! 



f (nS,o),- 








' (n;.„),- 




,(nS,i),- 


(n!a). 



-'-'-0,1 -'-'0,1 ^ -'-'1,1 -'-'-I 



0,0- 



2.34. Second auxiliary system. It is now time to come back to the functions G, H, L 
and M and to get rid of the auxiliary "Pi" functions. Unfortunately, we cannot invert 
directly the linear system (2.18), hence we must choose two specific square functions 
as principal unknowns, and the best, from a combinatorial point of view, is to choose 
and Dly. Remind that by (2.20), we have 0°^ = U^ q and = U\ ^. For clarity. 
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it will be useful to adopt the notational equivalences 
(2.35) e° = n°o and 6^ = Ul^. 

We may therefore quasi-inverse the linear system (2.18), obtaining that the four func- 
tions IIq q, nj ^, IIq ^ and 11° ^ may be expressed in terms of the functions G, H, L and 
M and in terms of the remaining two principal unknowns (2.35), which yields: 

= -G 

1° - 1 r -u 1 fii 

^° = M. 

yy 

Replacing now each of these four expressions in the compatibility conditions of the 
first auxiliary system (2.33), solving the four equations with respect to 6^, 6°, 6^ 
and Qy, we get after hygienic simplifications what we shall call the second auxiliary 
system, which is a complete system of first order partial derivatives in the remaining 
two principal unknowns 6° and 0^: 
(2.37) 









= □ 






nii = 


= □ 


(2.36) 


< 










■'■■'■0,1 " 


= □ 








= □ 





- Ly + 2M, + HM -II^ ^ 
^ 2 










e° = 


— 2 Gy + Hx -\- G L — — - 


-GQ^ 








ol = 


2 1 1 
3^3 2 


HL - 


2 2 






e° = 

y 


12 1 
- -Hy + -Lx + 2GM - - 
3^3 2 


HL - 


ii/e' + i 

2 2 


L0°H 


2 



We do not comment the intermediate computations, since they offer no new combina- 
torial discovery. 

2.38. Precise lexicographic rules. We group first order derivatives before zeroth order 
derivatives; in each group, we respect the lexicographic order of appearance given 
by the sequence G, H, L, M, 0°, 6^; we always put rational coefficient of every 
differential monomial in its left; consequently, we accept a minus sign just after an 
equality sign, as for instance in (2.36)i and in (2.37)2; for clarity, we prefer to write a 
complicated differential equation as = $, with on the left, instead of $ = 0, since 
$ may incoporate 10, 20 and up to 150 monomials, as will happen for instance in the 
next sections below. 

2.39. Compatibility conditions for the second auxiliary system. Clearly, the nec- 
essary and sufficient condition for the existence of solutions (0°, 0^) to the second 
auxiliary system (2.37) is that the two cross differentiations vanish: 



(2.40) 



Using (2.37), we shall see that we exactly obtain the two second order partial dif- 
ferential equations written in Theorem 1.2 (4). For completeness, we shall perform 
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completely the computation of the first compatibility condition (2.40) and leave the 
second as an (easy) exercise. 

First of all, inserting (2.37) and using the rule of Leibniz for the differentiation of 
a product, let us write the crude result, performing neither any simplification nor any 
reordering: 



fo=(e;.) -(ej), 



(2.41) { 



2 Gyy + H^y + Gy L + G Ly ^ H Hy — Gy — G Q] , + 0° 0?.+ 



+ 2 H^y 



'^L^^ -2G^M -2G M^ + ^H^L + ^H L^- 



1 



Next, replacing each first order derivative 6°, 6°, 6^, and 6^ occuring in (2.41) by its 
expression given in (2.37), we obtain (suffering a little) as a brute result, before any 
simplification (except that we put all second order derivatives in the beginning): 
(2.42) 

4 „ 2 







2G 



yy 



+ Gy L + G Ly - H Hy - GyQ^ + G Ly - 2 G - G H M+ 

2 V ; 2 ^ ^ 3 ^ 3 

+ 2GMe°- -iJLe°- ii/0°0i + -L{Qy + -{QyQ^- 
2 2 2 2 

-2G,^M -2GM,^ + ]^H^L + ]^H L^ + ]^H^ 0^+ 

+ \HL,~\HHy+GHM~- {HfL~ - [Hf 0^+ 
3 4 4 

+ -HLQ'^ + -HQ^Q^--L^Q^+GyL--H^L--G (Lf- 
4 4 2 ^2 2^^ 



G„0i-:,i7.0i- 



+ li72^ + lG'L0i--L (00)^ 
4 2 4 V ^ - 

-^GL0i + i/720^ + iG(0Y-l (0°)'0^-^L.0°+ 



4 



6 



2 

+ -i7y0°-G'M0° + -i/i:0° + -if0°0'--L (0°)2- 



3 
4 



4 



Now, we can simplify this brute expression by chasing every couple (or triple, or 
quadruple) of terms killing each other. After (patient) simplification and lexicographic 
ordering, we obtain the equation 







-2G. 



(2.43) 



yy 



in. 



xy 



2 

2 -^zx"t" 



2{GL)y -2G^M -AGM^ 



2 4 

— H Lt H Hy, 

3 3^' 



JOEL MERKER 



16 



which is exactly the first equation of (4) of Theorem 1.2. The treatment of the second 
one is totally similar. This completes the proof of the equivalence between (1) and (4) 
in Theorem 1.2. □ 



2.44. Interlude: about hand-computed formulas. In Section 4 below, when dealing 
with several dependent variables y^, . . . , y™, many simplifications of identities which 
are much more massive than (2.42) will occur several times. It is therefore welcome 
to explain how we manage to achieve such computations, without mistakes at the end 
and strictly by hand. One of the trick is to use colors, which, unfortunately, cannot 
be restituted in this printed document. Another trick is to underline and to number 
the terms which disappear together, by pair, by triple, by quadruple, etc. This trick is 
illustrated in the detailed identity (2.45) below, extracted from our manuscript, which 
is a copy of (2.42) together with the designation of all the terms which vanish together. 
Hence, we keep a written track of each intermediate step of every computation and of 
every simplification. Checking the correctness of a computation simply by reading is 
then the easiest way, both for the writer and for the reader, although of course, it takes 
time, anyway. 

On the contrary, when relying upon a digital computer, most intermediate steps are 
invisible; the chase of mistakes is by reading the program and by testing it on several 
instances. Alas, all the finest intuitions which may awake in the extreme inside of a 
long computation are essentially absent, the mind believing that the machine is stronger 
for such tasks. This last belief is in part true, in the case where straightforward known 
computations are concerned, and in part untrue, in the case where some new hidden 
mathematical reality is concerned. 

For us, the challenge is to control everything in a sea of signs. Computations are to 
be organized like a living giant coral tree, all part of which should be clearly visible in 
a transparent fluid of thought, and permanently subject to corrections. 

Indeed, it often happens that going through a problem involving massive formal 
computations, some disharmony or some incoherency is discovered. Then one has to 
inspect every living atom in the preceding branches of the growing coral tree of com- 
putations until some very tiny or ridiculous mistake is found. In addition to making 
easy the reading, a perfectly rigorous way of writing the formal identities which re- 
spects a large amount of virtual conventions facilitates to reorganize rapidly the coral 
tree after a mistake has been found. The accumulation of new virtual conventions, all 
of which we cannot speak, constitutes another coral meta-tree and another profound 
collection of trick. Finally, we use a blank fluid corrector to avoid copying to much. 



Extracted from our manuscript, here is the identity (2.42) with the underlining- 
numbering of all the vanishing terms (without the original colours) until we get the 
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final equation (2.43): 

4 2 

= — 2 Gyy + — Hxy — — Lxx-\- 

+ GyL+GLy- H Hy-Gy ^+ G L y -IGM^' GHM f^+ 

1 o n 1„ /c»ix2 ^ o- r" 2 



+ ^G (L? ^- GM 9° ^- - G (Q^? -i^HyO"^ + l,. 6° _+ 

^ © ^2 @ j3 

1 „ . 1 „ 



+ 2GM9° ^- -gLe° --/7G°e^ +-L(e°)2 +-(QOfQ^ - 



-2GxM -2GMx + \HxL+]-HLx + \Hxe^ + 
2 2 2 

+ \HLx-\HHy + GHM®- \ {Hf L - \ [Hf + 

— i\J — (m) 



(2.45) 



111 11 



i7Le° +-i70°0i --L,e° +GyL--HxL--G{Lf + 
4 ® i © 2 @ 2 2 @ 

4 © 2 © ^__© ® 2 Q 

2 © 4 @ 2 © 6 



3- ^© ^ i © i © ^ ® 

4 © 

As may be observed, the order in which we discover the terms which vanish is gov- 
erned by chance. After some terms are underlined, they are automatically disregarded 
by the eyes, which lightens the chasing of other terms to be simplified. To collect the 
remaining terms in order to obtain the final expression (2.43), our method is similar: 
we underline the terms which may be summed together. However, whereas we use the 
red pencil to underline the vanishing terms, we use the green pencil to underline the 
remaining terms. This small trick is to avoid as much as possible to copy several times 
some long formal expressions. Finally, we reorder everything lexicographically, so as 
to get the conclusion (2.43). In order to obtain the final equation (2.43) as efficiently 
as possible, we read the remaining terms, picking them directly in lexicographic order. 
If, by lack of luck, one or two terms are forgotten by the eyes and not written in the 
right place, we copy once more the very final result in the right order, or we use the 
blank corrector. 

Of course, such a refined methodology could seem to be essentially superfluous for 
such relatively accessible computations. However, when passing to several dependent 
variables, the current expressions will be approximatively five times more massive. We 
may really ascertain that a clever methodology of hand computations is helpful in this 
category. 
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§3. Systems of second order ordinary differential equations 
equivalent to free particles 

3.1. Combinatorics of the second order prolongation of a point transformation. In 

this section, we endeavour to explain how Lie's theorem and proof may be generalized 
to the case of several dependent variables. As in the statement of Theorem 1.7 (1), let 
us assume that the system yi^ = F^{x, y, y^), j = 1, . . . , m, is equivalent under an 
invertible point transformation (x, y) i— {X{x, y), Y{x, y)) to the free particle system 
^xx = 0, J = 1, . . . , m. By assumption, the Jacobian determinant 



(3.2) 



[x\y 



Xx Xyl 



Vm 

Si 



Xyr 
lyn 



•ym 

lyvn. 



does not vanish at the origin. As in the case m = 1, since the flat system Y^^ = 
is left unchanged by any affine transformation, we can (and we shall) assume that the 
transformation is tangent to the identity at the origin, so that the above Jacobian matrix 
equals the identity matrix at (x, y) = (0, 0), whence in a neighborhood of the origin it 
is close to the identity matrix, namely 



(3.3) 



X^, 



X, 



Inductive formulas for the computation how the differential equation in the {X, Y) co- 
ordinates is related to the differential equation in the (x, y) coordinates may be found 
in [IBK1989I , [IOL1995I ; the explicit formulas are not achieved in these references. Let 
us recall the inductive formulas, just on the computational level (differential-geometric 
conceptional background about graph transformations may be found in [I011986L 
Ch. 2). 

First of all, we seek how the Y^ := are explicitely related to the y^. It suffices 
to replace, in the identity 



(3.4) yI ■ Ix^dx + Y^ Xyl dy^ ) =Y^dX = dY^ = Y^ dx + Y^ Y^, dy' 



the differentials dy^ by yl dx and then to identify the coefficient of dx on both sides, 
which rapidly yields the formulas 



(3.5) 



Yi 



X 



^X + X]« = l Ux ^yl 

Xx + T.T=i y'xXy 



for j = 1, . . . , m. 



CHARACTERIZATION OF NEWTONIAN FREE PARTICLES 



19 



Next, we seek how the Y^j^ := = are related to the y^?, y''^^. It suffices to 
again replace each dyi by dx and each dy^ by y^^ dx in the identity 
(3.6) 



Y^^- (x^dx + f^^y^dy'] =Y^j,-dX = dYi 



1=1 



Before entering the precise combinatorics of the explicit expression of Y^X' 
observe that the last term of (3.6) simply writes D{Y^) dx, where D denotes the total 
differentiation operator (of order two) defined by 



Since dX = DX after replacing each dy'' by yl dx, it follows that we may compactly 
rewrite (3.6) as 

(3.8) Y^x DX-dx = D (Fi) ■ dx 
Consequently, the expressions of (obtained in (3.5)) and of Y^-^ are 

(3.9) yi=£r! and y> ^1 ^ ££XL^1^££L^. 

^ ^ ^ DX DX [DX]^ 

As, by assumption, the system y^^ = F^(x, y, y,j) transforms to the flat system = 
0, after erasing the denominator of (3.7), we come to the equations 

(3.10) = DDY^ ■ DX - DDX ■ DY^, 



for j = 1, . . . , m. However, this too simple and too compact expression of the system 
Vxx — F^{x, y, y^) is of no use and we must develope (patiently!) the explicit expres- 
sions of DDY^ , of DX, of DDX and of DY\ using the complete expression of D 
defined in (3.6). 

At this point, we would like to stress that it constitutes already a nontrivial compu- 
tational and combinatorial task to obtain a complete explicit formula for the system 
yi^ = F^{x, y, y.j) hidden in the compact form (3.10), which would be the generaliza- 
tion of the nice formula (2.9) involving modifications of the Jacobian determinant. For 
general m ^ 2, the complete proofs are postponed to Section 5 below. 

Since it would be intuitively unsatisfactory to provide directly the final simplified 
expression of the development of (3.10) in the general case m ^ 2, let us firstly de- 
scribe step by step how one may guess what is the generalization of (2.9). 
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For instance, in the case m = 2, by a direct and relatively short computation which 
consists in developing plainly (3.10), we obtain for j = 1,2: 
(3.11) 

0= -X.,1^+1^X.+ 



yl 



+ yl 



+ yl yl 



+ ylyl 



+ yl yl 



+ yl yl yl 



+ yl yl yl 



+ yl yl yl 



I 2 2 2 

+ yxyxyx 



-X.,, Yl + YL, - 2 Yl + 2 YL, 



-X,.. n + Xy. - 2 X.,. yi + 2Y^. X, 



+ 
+ 



-2 X^„i F{ + 2 X„i - y-?' + f{ 1 x^ 

J-y yi- ' xy^ y y y X ' y'-y^ 



+ 



-2 FJ, + 2 X,. - 2 X^y. Yl^ + 2 1^^, X,x - 

-2x,vn^' + 2i;;,„,x, 

-2 X^.„2 Y\ +2Y\ X„2 - X^2^2 !?■ + 2 X^. 

xy y^ I ^.^2 y y y X ' y'y^ 

— X..1..1 Y'^i -\- Y\ -1 X,,i -|- 
y y y y y y 

-Xy.y. Yl, + Yl^^, Xy. > 2 X,.,2 F^, + 2 rj,^, X,. 
-X,2,2 + Yl^^, Xy. - 2 X,.,2 Yl, + 2 rj,^, X,2 



+ 



"X,,2,,2 y^"'o -|- yl.^o X,, 



+ 



yxx 



a XX 



-Xy. Fi + FJ, X. + ^/^ <! -X,. + Yl X, 



+ 



-X,2 Yi + FJ, X. + y^ { -X,2 Yl, + FJ, X,i 



Unfortunately, above two equations are not solved with respect to yl^ and to y. 
Consequently, if we abbreviate them as a linear system of the form 



2 

XX' 



(3.12) 



= A' + yl-Bl + yl-Bl, 
= A' + yl-Bl + yl.Bl 



we have to solve for yl^ and for y^^ by means of the classical rule of Cramer. Here, it 
is rather quick to check manually that the determinant of this system has the following 
nice expression: 



(3.13) 



Bl Bl 
Bl Bl 



= Ai 
= Ai 



X 



\y^\y'')-{X^, + ylXy.+ylXy2} 



{x\y'\y^)-DX. 



However, the complete solving for yl^ and for y^^ requires some more time. After a 
direct and rather long hand computation (or alternately, using Maple or Mathematica) 
one obtains formulas involving hidden 3x3 determinants, which have to be guessed 
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(3.14) 



Vxx 



yl 



+ yi 



X 



Xyl 



Xy2 

yl 

y^ 

y\ 

y^ 



Xx 

yl 

X 

Y3 



Xxx 

yl 

XX 



yl 



+ 





^xj/-'- 


Xy2 


yl 


yl 

xy^ 

Y'\ 

xy 


yl 

iy2 

y\ 

y^ 






X 2 


Xy2 


yl 

X 


Y\ 

xy' 

Y'\ 

xy^ 


yl 

^y' 

y2 


y2 

X 



X 

yl 

XX 

y2 

XX 



+ 



Xyl 

yl 

yl 

y^ 



Xy2 

yl 

y| 



+ 



+ ylyl 



I 2 2 
+ ?/x 1/x 



Xylyl 




yl 
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yiyi 

loo 

y^y^ 



Xy2 

Y\ 

y^ 

y2 

^y2 
Xy2 

Y\ 

/2 

iy2 
Xy2 

yl 

-f y2 

y| 



- 2 



+ 



^xy^ 




Xj,2 


Y\ 


y\ 


Y\ 


y^ 


yt 

xj/i 


y2 


y2 





Xyl 


Xy2 


Y\ 


Y\ 

y^ 


y^2 




y^ 


y^2 


y2 


y2 








+ yl yl yl ■ 



+ yl yl yl 



Xylyl 

yl 

Y\ , 

yl 

y y 
y2 

y2y2 



Xyl 

yl 

^y^ 

Y\ 

y^ 

Xyl 

Y\ 
y^ 
y2 

^y^ 



Xy2 

yl 

ly2 

y\ 

y 

Xy2 

Y\ 
y^ 

Yl 



+ yl yl yl 



Xyly2 


Xyl 


Xy2 


Y\ 2 

yLyA 


yl 

^y^ 


yl 

iy2 




Y\ 2 

y y 


Y\ 

y 


y\ 


y 



+ 



This formula and the next have been checked by Sylvain Neut and Michel Petitot with 
the help of Maple. We notice that the size is not negligible, but fortunately, there 
appears some combinatorics, much more visible than in (3.11). The second equation 
that we obtain, namely for y^^,, is as follows: 



(3.15) 



+ yl 



+ yl 



Xyl 

Y\ 
y^ 
y2 

^y^ 



yl 

X 

y2 



y 



Xy2 

Y\ 

y 

Yl 



Xx Xyl Xx 

yl 

x 

y2 



Xx Xyl Xxyl 

yl yl 

j/i xj/i 

y2 y2 

j/l xj/l 

Xx Xyl Xxy2 

yl yl yl 

y^i y^^ 

s/ xy 



yl yl 

-"j/i -"xx 

y| y2 

y2 XX 



+ 



+ 



^XX 




Xy2 


yl 

XX 


yl 


yl 

-f „2 


y2 

XX 


y^ 

y^ 


y^ 

y^ 



+ 
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+ ylyl 



+ ylyl 



I 2 2 

+ yxyx 



y2 



Xyl 



^yl ^ylyl 

y y y 



Xx Xyl 

Yl Y^\ 
Y^ Y\ Y\ 2 



Xyly2 

Y\ 2 

yLyJ. 



y' y'i y| 2 

x y2y2 



y2y2 




+ yl yl yl 



+ yl yl yl 



Xylyl 

yl 

y2 
^yW 

y2y2 

Yl., 
yy 



Xyl 

yl 

y2 

Xyl 

yl 

^y^ 

Y\ 

y^ 



Xy2 

yl 

^y2 

y2 

yl 

-f „2 

y| 



+ y 



Xyly2 

Y\ 2 

yly2 

Y\ . 

yly2 



Xyl 

yl 

y2 



Xj^2 

yl 

ry2 
Yl 



+ 



Importantly, the obtained formulas seem to be analogous to the formula (2.9), since we 
observe that the coefficients of the degree three polynomial in the are modifications 
of the Jacobian determinant A{x\y^\y'^). 

To describe the underlying combinatorics, let us observe that there exist exactly six 
possible distinct second order derivatives: xx, xy^, xy^, y^y^, y^y^ and y^y^. There 
are also exactly three columns in the Jacobian determinant (3.2). By replacing each 
of the three columns of first order derivatives by a column of second order detivatives 
(leaving X, yi and Y'^ unchanged), we may build exactly eighteen new determinants 



(3.16) 



' Aixx\y'\y^) 


A{x\xx 


y') 


A{x 


y^\xx) 


A{xy'\y'\y') 


A{x\xy'^ 


y') 


A{x 


y^\xy^) 


Aixy'\y'\y') 


A{x\xy'^ 


y') 


A{x 


y^\xy'^) 




A{x\yY 


y') 


A{x 


yW 




A{x\yY 


y') 


A{x 


yW 


. HyW\y') 


A{x\yY 


y') 


A{x 


yW 



where for instance 



(3.17) 





Xyly2 


Xyl 


Xy2 




A{yW\y') := 


Y\ 2 

y y 


Y\ 

y^ 


Y\ 
y 


and 




Yl 2 

y y 


Yl 

y^ 


Yl 

y^ 





A{x\y \xy ) : 



X-x -Xyl Xxy2 

yl yl Y~ 

^x -fyl J^xy^ 

Yl Yl Y\ 

X xy^ 



Hence, using the A-notation, we may rewrite the two equation (3.14) and (3.15) under 
a more compact form; after division by the Jacobian determinant A{x\y^\y'^), the first 
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equation becomes: 
(3.18) 



+ yl 



+ yi 



+ ylyl 



+ ylylyl 



+ yl yl yl 



^W\y') 

A(x|x?/^|?/^ 



A(x|xy^|y^) A(a^aVb^l 



A(x|i/i|i/2) 

A(x||/V|y^) , 
A(x|yi|2/2) ^ 

A(yV|y^|l/^) 
A(a;|2/%2) 

A(yV|y^|y^) 
A(x|2/%2) 



A(x|y%2) 

A(V|yi|y^) 
A(x|yi|y2) 



l/x l/x 1/x 



2 2 



A(x|yi|y2) 

A(x|2/V|y') 



A(x|yi|y2) 
A(yybV) , , 



Similarly, the second equation takes the form: 
(3.19) 

2 , A{x\y^\xx) 1 f A(x|?/Vz/^y 



= y^. + 



+ yl 



A(x|yi|y2) 



A(x|?/^|x?/^) A(xx|y"'"|y 



A(x|i/i|y2) 



A(x|i/ib2) A(x|yi|y2) 



yl yl 



I 2 2 

+ yxyx 



+ ylylyl 



I 2 2 2 

+ yxyxyx 



A(x|i/%y^ 



A(xyVV) 
A(x|2/i||/2) 

A(VbV)l 
A(a;|yi|y2) A{x\y^y') j 

A{yY\y'\y')\ , ^ , , 

yxyxy. 



A(x|y%2) 
A(a;|2/2|2/V) 



A{x\y^\y^) 
A{x\y^\y^) 



A{x\y^\y^y^\ 
A(x|y%2) 



-2 



A(yVbV)l 



Since the formulas are still of a consequent size, analogously to what was achieved 
in Section 2, we shall introduce a new family of square functions as follows. We 
first index the coordinates (x, y^, . . . , y™) as y^, . . . , y*"), namely we introduce 
the notational equivalence 



(3.20) 



X 



which will be very convenient in the sequel, especially in order to write general for- 



mulas. With this convention at hand, our eighteen square functions D\ , defined for 



y '-y ■ 
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^ ji, j2, ki ^ 2 axe defined by 
(3.21) 



□0 

— XX 


A{xx\y^\y'^) 
■ A(x|yi|y2) 

\ It? It? J 


□° 1 


A(xyi|yi|y2) 
■ A(x|?/i|y2) 

\ It? \iJ J 


□° 2 


A(xy2|yl|y2) 

■ A(x|yi|y2) 

\ It? 1 i? / 




A(2/iyi|yi|y2) 
A(x|yi|2/2) 

\ It? \iy / 


1— lO 

□ 1 2 


_ A(yV|y%') 

Afxlyilu^) 

\ It? It? / 




A(y2y2|yl|y2) 

A(x|i/My2) 

V It? 1 £? / 


' — 'xx 


A{x\xx\y'^) 


^ly^ 


_ A(x|xy^|y2) 

AfxIuMy^) 

V 1 y 1 y / 


□ ^ 2 


_ A(x|xy2|?/2) 
A(x|?7i|?y2) 


< 

^y'y' 


A(x|yi2/i|y2) 
• A(x|yi|y2) 


□Jlj^2 


A(x|yiy2|y2) 
• A(x|yi|y2) 


□J2j^2 


A(x|y2y2|y2^ 

• A(x|yi|y2) 




A(j; y-'- xx) 


1 


_ A{x\y^\xy^) 


2 


_ A(x|yV2/^) 


— XX 


A(x|yi|y2) 




A(x|2/i|2/2) 




A(x|yi|y2) 


^y'y' 


A(x|y%iyi) 
■ A(x|yi|y2) 




A(x|yi|yiy2) 
■ A(x|yi|y2) 


n2 


A(x|yi|y2y2) 
■ A(x|yi|y2) 



Obviously, the square functions are symmetric with respect to the lower indices: 
^yHyi2 ~ ^yhyh- Hcrc, thc uppcr index designates the column upon which the sec- 
ond order derivative appears, itself being encoded by the two lower indices. Even 
if this is hidden in the notation, we shall remember that the square functions are 
explicit rational expressions in terms of the second order jet of the transformation 
(x, y) I— s> y), Y{x, y)). At this point, we may summarize what we have estab- 

lished so far. 

Lemma 3.22. The system of two second order ordinary differential equations y].^ = 
F^{x, y, yx) and y1^ = F^{x, y, y^) is equivalent, under a point transformation, to the 
flat system Yxx = and Y^x = i/ and only if there exist three local K-analytic 
functions X{x, y), y) and Y'^{x, y) such that it may be written under the form 

(3.23) 

' = 2/L + nL + 2/^ (2 n^,! - + 2/^ (2 □^,2) + 

+ ylyl- [Ul^y, - 2 Uly,) + ylyl ■ (2 □J,^2 - 2 □"^2) +ylyl- (□^2,2) + 

+ ylylyl ■ {-^W) + vlvlvl ■ {-'^^W) + vlvlvl ■ (-^S^yO > 
' = + nL + z/x ■ (2 + yl ■ (2 uly, - + 

+ ylyl- {Ul^y,) + ylyl ■ (2 □^,^^2 - 2 Uly,) +ylyl- (□^2,2 - 2 □"^2) + 

+ yl yl yl ■ {-^yw )+ylylyl- (-2 nSi.O + yl yl yl ■ (-^S^yO • 

3.24. Second Lie prolongation of a vector field. At this point, instead of proceeding 
further with the case m = 2, it is now time to pass to the general case m ^ 2. First of 
all, we would like to remind from nGM2003ll the complete explicit expression of the 
point prolongation to the second order jet space of a general vector field of the form 
L = X # + V™ , #7: it is a vector field of the form 




where the coefficients R{ and are polynomials in the jet space variables having as 
coefficients certain specific linear combinations of first and second order derivatives of 
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X and of the F^ : 



(3.26) 



m 

^xx ~^ ^ ^ Vx 

m 



m m 



«1 = 1 l2=l 



[-61 Xyl^] 



2Y^ 



6u + 



m m 



yhyh 



xy'- 



SIX. 



h xy' 



+ E E E y'^' y'^ ■ Hi ^y^^y^-^] + 

h=l l2 = l l3 = l 
m 

+ E^^-[^.-2^^^^ 

ll=l 

m m 

+ EE^'^^-K^. 



ll = l l2 = l 

However, since the notations in nGM2003ll are different and since the general case of 
n ^ 1 independent variables and m ^ 1 dependent variables is considered there, it 
is certainly easier to reconstiture formulas (3.26) directly by means of the inductive 
formulas described in I1011986L [|BKl989l). 

Analogously to the observation made in Section 2, we guess that there exists a formal 
correspondence between the terms of R2 not involving y^^ and the explicit form of the 
equation yi^ = F^{x, y, y^) equivalent to Y^-^ = 0. In the case m = 2, we claim that 
this formal correspondence also holds true. Indeed, it suffices to write formula (3.26) 
for R2 modulo the y^^,, which yields two expressions in total analogy with the two 
explicit polynomials appearing in the right-hand side of (3.23): 
(3.27) 



( (mod yiJ = + 1/M2 Y^y. -Xxx]+yl-{2 Y^y, ]+ylyl- - 2 X, 



yl 
ylvl 



'-xy'- 



2 Xxy2 J" 



ylyl ■ {^j/v} 



+ yl yl yl ■ {-^y^y^} + yl yl yl ■ {-2 Xyiy2} + yl yl yl ■ {-Xy2y2 



R2 



mod y] 



XX I 



Yl 



yx ■ 



{2 Yly^ } + Z/M2 Yly2 - + yl yl ■ [Y^.y. } 



+ ylvl ■ {2 Ylry2 - 2 X^y. ]+ylyl- {Yy%2 - 2 X,y2 } + 

+ yl yl yl ■ {-^y^y^ } + vlvlvl ■ {-'^ ^yw } + ylylyl- {-^; 



Except for inductive inspiration {see the formulation of Lemma 3.32 below), this ob- 
servation will not be used further. At this stage, it helps at least to maintain a strong 
intuitive control of the correctness of the underlying combinatorics. 

3.28. System equivalent to the flat system. By induction, we therefore guess that 
the analogy holds for general m ^ 2, namely we guess the following combinatorics, 
which requires some preliminaries. 

As in the beginning of §3.1, let x G K, let y = (y^, . . . , y™) G K™, let (x, y) t-^ 
y), Y(x, y)) be a local K-analytic transformation defined in a neighborhood of 
the origin in K*""*"^ and assume that the system yl^ = F^(x, y,yx), j = 1, . . . , m, is 
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equivalent to the flat system Y^j^ = 0, j = 1, . . . , m. By assumption, the Jacobian 
matrix of the equivalence equals the identity matrix at the origin. Remind that we 
identify x with y^. For all ki, /i, /2 = 0, . . . , m, we define a modification 

(3.29) A(x|...|Vi/^^|...|y") 

of the Jacobian determinant as follows. We replace the /ci-th column of the determi- 
nant (3.2), which consists of first order derivatives -yk^ , by a column which consists of 
second order derivatives -yhyh ■ In (3.29), the notation 1'^^ designates the fci-th column, 
the first one being labelled by A;i = and the last one by ki = m. With this notation 
at hand, we may define the square functions 



(3.30) 



r-|fci 

^yhyh 



A{x 




k 








A( 


X 




kiyki 1 


• --y 





which are rational expressions in the second order jet of the transformation (a;, y) ^ 
{X(x,y),Y{x,y)). As before, the denominator is the Jacobian determinant of the 
change of coordinates. 

Since, according to (3.26), the expression of R2 (mod y^^) is 
(3.31) 



R^2 (mod yl 



h=i 



+ 



yhyh 



cyh 



h = l l2 = l 
m m m 

h = l «2 = 1 ^3 = 1 



y'.'y'^y'^-[-5lXy.,y.,] 



+ 



and since, in the cases m = 1 and m = 2, we have already observed strong analogies 
between (3.31) and the complete explicit expression of the system yi^ = F^(x, y, y.^) 
equivalent to the flat system Y]^-^ = 0, we guess that the following lemma is formally 
true. 

Lemma 3.32. The system y^^ = F^{x, y, y^), j = 1, . . . , m, is equivalent to the flat 
system Y^-^ = 0, j = 1, . . . , m, if and only if there exist local M.- analytic functions 
y) and Y^{x, y), j = 1, • • • , m, such that it may be written under the speciflc 

form 



(3.33) 







Vxx 



m. 



f 

XX 



+ 



m m 



Zl=l l2 = l 



+ 



+ yiJ2J2y- 

Zl=l l2=l 



y'x' 



^yhyh 



The complete proof of this lemma involves only linear algebra considerations, al- 
though with rather massive terms. This makes it rather lengthy. Consequently, we 
postpone it to the final Section 5 below. 
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3.34. First auxiliary system. Clearly, if we set 



(3.35) 



K 

Li I 



xy' 



M, 



yhyk 



we immediately see that the first condition of Theorem 1 .7 holds true. Moreover, we 
claim that there are m + 1 more square functions than functions G^, H^^, Lj_^ and 
Mi^^i^. Indeed, taking account of the symmetries, we enumerate: 

#{□1} = 1, 



(3.36) 

whereas 
(3.37) 



m 



im?{m + 1) 



H'i^lhyh) 



m(m 



m?{m + 1) 



mm 



m 



1) 



Similarly as in Section 2, for j, Zi, I2 = 0, 1, 



m 



, let us introduce functions of 
{x,y^, . . . , y™), symmetric with respect to the lower indices, and let us seek necessary 
and sufficient conditions in order that there exist solutions {X, Y) to the first auxiliary 
system defined precisely by: 

^xx — -^^0,0' ^xyh ~ -^^0,/i' '-'v'lVa ~ "^^^1,^2' 



(3.38) 



U,U' xv'l U,(l' yny'2 11,12 



xy'i 



-yhyh 
yhyh 



3.39. Compatibility conditions for the first auxiliary system. As in Section 2, the 
compatibility conditions for this system will simply be obtained by computing the 
cross differentiations. The following statement generalizes Lemma 2.31 and also pro- 
vides a proof of it, in the case m = 1. 

Lemma 3.40. For all j, /i, I2, /a = 0, 1, . . . , m, we have the cross differentiation rela- 
tions 



(3.41) 



yhyh 



y^y 



k=0 k=0 

Proof. To begin with, as a preliminary, let us generalize the Pliicker identity (2.28). Let 
Ci, C2, . . . , Cm, D, E he (m + 2) column vectors in K"* and introduce the following 
notation for the m x (m + 2) matrix consisting of these vectors: 



(3.42) 



[Ci|C2|---|C^|/^|^] 



Extracting columns from this matrix, we shall construct m x m determinants which 
are modification of the following "fundamental" determinant 



(3.43) 



||C*l| ■ ■ ■ |C*m|| = He'll ■ ■ ■ ^"'"C'j-^l ■ ■ ■ P^Cjjl ■ ■ ■ \Cri 



Here and in the sequel, we use a double vertical line in the beginning and in the end 
to denote a determinant. Also, we emphasize two distinct columns, the ji-th and the 
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j2-th, where j2 > ji, since we will modify them. For instance in this matrix, let us 
replace these two columns by the column D and by the column E, which yields the 
determinant 

(3.44) \\Ci\---\^W\---\^^E\---\Cm\\. 

In this notation, one should understand that only the ji-th and the j2-th columns are 
distinct from the columns of the fundamental m x m determinant (3.43). With this 
notation at hand, we can now formulate and prove a preliminary lemma that will be 
useful later. 

Lemma 3.45. The following quadratic identity between determinants holds true: 
(3.46) 





p'^l ■ ■ 


rE\---\Cn 






hn. 1 . . . 1 . . . 


\Cn\\ - 


= He'll 


■ • ■ 


. . . par'. 1 . . . 


|Cn|| ■ 




■■■\^-'Cj,\---\^^E\ 


■■■\Cn 


- He'll 


■ ■ ■ \''E\ 


|i2^j2 1 . . . 


|Cn|| ■ 




■■■rc,,\---rD\ 


■ ■ ■ \Cn 



Proof. After some permutations of columns, this identity amounts to 

||Cl| ■ ■ ■ |C*m-2|-D|-E'|| ■ ■ ■ ■ |Cm„2|Cm-l|C*m|| = 

(3.47) { = IICll ■ ■ ■ |C'm-2|-D|Cm|| ■ ||Ci| ■ ■ ■ |Cm_2|C'm-l|-£^|| — 

— ||Ci| ■ ■ ■ |Cm_2|-£'|C'?7i|| ■ He'll ■ ■ ■ |C'm-2|C'm-l|-D|| • 

To establish this identity, we introduce some notation. If A and B are vertical vectors 
in and if ii,i2 = I, . ■ ■ ,m with ii < i2, we denote 



(3.48) K,M\B) 



Ai^ Bi-^ 
A ■ R- 

^«2 -'-^«2 



If ||y4i|y42|v43| ■ ■ ■ \Arr\ is a m X m determinant, and if ii,i2 = 1, . . . ,m with ii < 12, 
we denote by M^'^^A^^] ■ ■ ■ \ Am) the (m — 2) x (m — 2) determinant obtained from 
the matrix [Aal • ■ ■ 1^1^] by erasing the ii-th line and the ?2-th line. Without proof, we 
recall an elementary classical formula 

( \\A^\A2\A^\---\Am\\ = 
^^•"^^^ I = E i-^y'^''''K.MM2)-M-T^\A,\...\AJ, 

which may be established by developing the determinant ||yli|A2|A3| ■ ■ ■ \Am\\ with 
respect to its first column, and then re-developing all the obtained (m — 1) x (m — 1) 
determinants with respect to their first columns. To establish (3.47), we start with an 
equivalent version of the identity (2.28): 



(3.50) 



K..iD\E) . AlJC,\C2) = ^,,,{D\C2) ■ Al,^{C,\E)- 



where 1 ^ zi < Z2 ^ and 1 ^ ^3 < ^4 ^ m. Multiplying by i)*i+*2+»3+«4 2^ 
multiplying by MJ^J^^Cs] ■ ■ ■ \CJ, and multiplying by MJ^J^'iC^] ■ ■ ■ \CJ applying 
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the double summation Ei^i,<i,<i,n El<i^,<u^m' we get 
(3.51) 



M™7,^(C3|---|C™)-M™7^(C3 



\Cri 



E E 

-A^^,,^(E|C2) ■ Km(Ci\D)] ■ M-7^2(C3| ■ ■ ■ \CJ ■ M™-/(C3| ■ ■ ■ \CJ. 

Thanks to the relation (3.49), this last identity coincides exactly with the desired iden- 
tity (3.47). The proof is complete. □ 

We can now establish Lemma 3.40. As a preliminary observation, by the Leibniz 
rule for the differentiation of a determinant, we must differentiate every column: 

[A {y^'y'^'l- ■ ■ 1 2/''" 2/^'")]^, = A {y'y^'y'''\ ■ ■ ■ W^y''"') + ■■■ + 

+ A(y'^/i|---|i/V'"/'")- 
Using also the rule for the differentiation of a quotient, we may endeavour to compute 
the cross differentiations ( 



(3.52) 



of the left-hand side of (3.41). 



This will generalize (2.25). Sometimes in the computation, we shall abbreviate the 
Jacobian determinant A(y°| ■ • • |y"') using the shorter notation A; as before, a product 
between two elements of K will often be denoted by the sign for clarity. Here is 
the computation: 
(3.53) 

yhyh 



9 
















■ ■ ■ ' 





8 






y^^ 








{ A(yc 


■ ■ ■ ' 


r) 



A(yV^i 
+ A(y° 



\'y'W'\-- 
h''y''y''\ 



ly™) ■ A + ■ ■ ■ + 



A (/I 
A (/I 



W'y"^) ■ A- 

b") ■ [A [y^^ 



A (/I 



+ 



A(yV 
+ A(y° 



\lyhyhyh\ 



|y'")-A + --- + 
■■■r)-A^+- 



A(y°| 
A(y°| 



bS'") ■ A- 

ly'")- [A (yVl 



-A [y 



2/™)+--- + 



Crucially, we observe that all the determinants involving a third order derivative upon 
one of their columns kill each other and disappear: we have underlined them with @ 
appended. However, it still remains plenty of determinants involving a second order 
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derivative upon two different columns. We must transform all of them and express 
them in terms of determinants involving a second order derivative upon only one col- 
umn. To this aim, as an application of our preliminary Lemma 3.45, we have the 
following relations, valid for ji, j2, h^hjsji = 0, ■ ■ ■ and ji < j^: 
(3.54) 

r A (/I ■ ■ ■ p-^yS'^l ■ ■ ■ rvVl ■ ■ ■ ly'") ■ A (/I ■ ■ ■ ry^^l ■ ■ ■ ■ ■ ■ b'") = 





. . \i'iyh\^ . 






..\nyn\. 




■ ■ 


\...\nyhyh\. 






A [y'\ ■ 


.\nyh\. 







With these formulas, we may transform the lines number 3, 4, 5 and 8, 9, 10 of (3.53). 
Also, we observe that the lines 6, 7 and 11, 12 of (3.53) involve determinants having 
a single second order derivative. Taking account of the |^ factor, we deduce that the 
lines 6, 7 and 11, 12 of (3.53) may already be expressed as sums of square functions. 
Achieving all these transformations, we may rewrite (3.53) as follows 
(3.55) 

'a 



[A]^ 



A(2/V|---|V|---|2/") ■ A(yr ■ 
-A(y'^y'^|---|V|---|l/'")-A(y« 
+ ■ ■ ■ + 

+ A(y°| 
-A(y°| 



\,yhyh\ 



\3y^yh 



\Jyhyh\ 
\Jyhym.\ 



|y'")-A(y°| 



\'y'\---\yS'^)- 
\'y'\---\y'W') 



k=0 



Iv"") ■ A {y' 



\Jyhyh\ 



V|---|i/™)-A(?/°|---lW^I---b"^) + 



A(yVl---|V 
-A(y'i/^|--- 
+ ■■■ + 

+ A {y'\ ■ ■ ■ \^y^^y'-^\- ■ ■ \y^) ■ A ■ ■ ■ \^y^ ■ ■ ■ 

- A (/I ■ ■ ■ Py'^y™! ■ ■ ■ \y^) ■ A (/I ■ ■ ■ Py^l ■ ■ ■ \y'^y'^) 



'2 j/fe • 



k=0 



Notice that two pairs of "cdots" terms + • — h appearing in the lines 3, 4 and 8, 9 
of (3.53) are replaced by a single "cdots" term + ■ ■ ■ + in the lines 4 and 9 of (3.55). 
Importantly, we point that in the "middle" of the two "cdots" terms + ■ — h appear- 
ing in the lines 4 and 10 of (3.55) just above, there are two terms which do not occur: 
they simply correspond to the two underlined terms having (a) appended appearing in 
the lines 4 and 9 (3.53). 
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Now, taking account of the factor we can re-express all the terms of (3.55) as 
sums of square functions: 



(3.56) 



fc=0; k^j fc=0; k^j 

m 

^^yhyh ^yhyk- 



k=0 
m 



^%y'=^ihyh+ X] ^lhyh^ll2yk + 

k=0; kj^j fc=0; k^j 

m 



k=0 



Finally, we observe that in the two pairs of sums having k ^ j appearing in the lines 2 
and 4 just above, we can include the term A; = j in each pair, because these two terms 
are immediately killed inside the corresponding pair. In conclusion, after a final obvi- 
ous killing of four (among six) complete sums in this modification of (3.56), we obtain 
the desired formula (3.41), with two sums. This completes the proof of Lemma 3.40 
and also at the same occasion, the proof of Lemma 2.31. □ 

3.57. Compatibility conditions for the first auxiliary system. According to the (ap- 
proximate) identities (3.3), taking account of the explicit definitions (3.30) of the 
square functions, we have 



(3.58) 



\y = 

' — 'xx 

— XX 



^ XX 1 



^xyh 



xyh 



^yh 
V 1 



h yhyh' 



Consequently, the first auxiliary system (3.38) looks approximatively like a complete 
second order system of partial differential equations in the (m + 1) independent vari- 
ables (x, y) and in the (m + 1) dependent variables {X, Y). By means of elementary 
algebraic operations, one may transform this system in a true second order complete 
system, solved with respect to the top order derivatives, namely of the form 



(3.59) 



yj 

XX 



^^0,0' 
^^0,0' 



Xxyh 

yi 

xyh 



- A° 
A-'' 



Xyh 



y^y ■ 



Yyh 



y'ly". 



A^' 



where the A'^-Ij^ are local K-analytic functions of {x,y^\ X,Y^ , X^, Xyij,Yj ,Y^,J. 
For such a system, the compatibility conditions [which are necessary and sufficient for 
the existence of a solution {X, Y)] follow by obvious cross differentiation. Coming 
back to the system 3.38, these compatibility conditions amount to the quadratic-like 
compatibility conditions expressed in Lemma 3.40. In conclusion, we have proved the 
following intermediate statement. 



Proposition 3.60. There exist functions X, Y^ solving the first auxiliary system (3.38) 
of nonlinear second order partial differential equations if and only if the right-hand 
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side functions i^{x, y) satisfy the quadratic compatibility conditions 

(3-61) -^-^ = -Y. Ki. ■ + E Kis ■ K. 



Qyh Qyh 

^ ^ fc=0 fc=0 



forjMMM = 0, 1, . . . ,m. 

3.62. Principal unknowns. As there are (m + 1) more square (or Pi) functions than 
the functions G\ Hl^, and M;^ ^2 defined by (3.35), we cannot invert directly 
the linear system (3.35) (which is of maximal rank). Hence we must choose (m + 1) 
specific square functions, calling them principal unknowns, and similarly as in §2.34, 
the best choice is to choose and D^^,, for j = 1, . . . , m. For clarity, it will be 
useful to adopt the notational equivalences 

(3.63) e° = n°o and 6^' = I^j. 

Then we may quasi-inverse the system (3.35), which yields : 
(3.64) 

■'■■'■0,0 ~ '-'xx — 5 

TT^ - n-' - -- + - o° 

m , = . = -u , + i 5/ + - 5/ L\\ + - 5/ + - 5/ e's 

^^0,/i — ^xyh — 2 h,li 2 

Before replacing these new expressions of the functions IIq q, Hq^^, 11;^^ H^^ and 
Ilg into the compatibility conditions (3.61), it is necessary to expound first (3.61), 
taking account of the original splitting of the indices in the two sets {0} and 
{1,2,..., m}. This yields six families of compatibility conditions, totally equivalent 
to the compact identities (3.61): 
(3.65) 

f mm 

(n^o.o) , , - (n^o,0. = -"0,0 Ko - E <o ni,, + K,^ n^o,o + E <h n^, 

^ ^ k=l k=l 

m m 

{Km)^ - {%,o) = -<i. Ko - E + <o K,o + E <o ni„ 

^ ^ k=l k=l 

m m 

- (Kil,. = -<i. Ko - E <h Kk + Kis Ko + E <h "L., 



y'3 \ ---/J/' 



k=l k=l 



^ k=i ^ fc=i 



m m 



(n?„J, - K,o),.. = -Km Ko - E nf,,., K,k + n?,on° + E KA, 

k=l k=l 
m m 

iKh)yh ~ iK,h)yi2 = ~K,i2 Kfi ~ E Km K,k + K,h K,o + E K,h K,k- 



k=l k=l 
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3.66. Convention about sums. Up to the end of Section 4, we shall abbreviate any 
sum J2k=i J2k Sp- Such sums will appear very frequently. For all 

other sums, we shall precisely write down the domain of variation of the summation 
index. 



3.67. Continuation. Thus, we have to replace (3.64) in the six identities (3.65). 
Firstly, let us expose all the intermediate steps in dealing with the first identity (3.65)i. 
Replacing plainly (3.64) in (3.65)i, we get: 
(3.68) 



'■0,0 



y 



2 'i 2 '1 



Ci) k k Cb) 

Eliminating the underlined vanishing terms with the letters a, h, c and d appended, 
multiplying by —2 and reorganizing the identity so as to put the term 5j_^ 9° solely in 
the left-hand side, we obtain the relation 



(3.69) 



k k 

- \ E Hi - 51 + 1 51 e° 



3.70. Conventions for simplifications of formal expressions. Before proceeding fur- 
ther, let us explain how we will organize the computations with the formal expressions 
we shall encounter until the end of Section 4. Our main goal is to devise a method- 
ology of writing formal computations which enables to check every computation visu- 
ally, without being forced to rebuild any intermediate step. In fact, it would be 
unsatisfactory to just claim that Theorem 1.7 (3) follows by hidden massive formal 
computations, so that we have to guide the rigorous and demanding reader until the 
very extremal branches of our coral tree of formal computations. 

As an example, suppose that we have to simplify the equation Q = + By + A — 
B + 2C~\D-lA + \D + E + B- 2C.lni\\Q beginning, the terms A,^. and By 
are differentiated once and they do not simplify with other terms. To distinguish them. 
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we underline them plainly and we copy the nine remaining terms afterwards: 

r = A^+ By + 

Here, each remaining term is also underlined, with a number or with a letter appended. 
For reasons of typographical readability, we never underline the sign, + or — of each 
term; however, it should be understood that every term always includes its (not under- 
lined) sign. Until the end of Section 4, we shall use the roman alphabetic letters a, h, 
c, etc. inside an octagon Q to exhibit the vanishing terms. As readily checked by the 
eyes, we indeed have —B_qj\- B_q= and 2 C ^^— 2 C (^= 0. Also, until the end of 
Section 4, we shall use the numbers 1, 2, 3, etc. inside a square □ to exhibit the remain- 
ing terms, collected in a certain order. The numbers have the following signification: 
after the simplifications, the equation (3.71) may be written 

r = + By + 

O-il) { 11 

+ - A- - D + E. 
\ 3 6 

Here, the plainly underlined terms + By do not count in the numbering (their num- 
ber is zero, for instance) and the first term of the second line | A correspond to the 
addition of all terms \T\ in (3.69). Analogously, the second term — | D correspond to 
the addition of all terms [2] in (3.69). Again, this guiding facilitates the checking of the 
correctness of the computation, using simply the eyes. No hidden delicate computa- 
tional step is "left to the reader" for the convenience of the writer. 

This principle will be constantly used until the end of Section 4; it has been system- 
atically used in IIM2004II and it could be applied in various other contexts. Again, the 
advantage is that it enables to check the correctness of all the formal computations just 
by reading, without having to write anything more. This is also useful for the author. 

3.73. Choice of an ordering. Until the end of Section 4, we shall have to deal with 
terms G, H, L, M, 6 together with indices and partial derivatives up to order two. In 
order to organize the formal expressions in a way which provides an easier decipher- 
ing, it is convenient to introduce an order between these differential monomials. In a 
symbolic index-free notation, we choose: 

(3.74) G < H < L < M <e. 

It follows for instance that G < GH < GL < HHL < HLMQ. Also, if a sum 
appears, we choose: GM < ^ GM. 

Here, we have only considered terms of order zero, without partial differentiation. 
The first order partial differentiations are and {■)y, again in symbolic notation, 
dropping the indices. We choose: 

(3.75) G^ < Gy < < Hy < < Ly < < My < < Qy < G < ■ ■ ■ . 
For second order derivatives, we choose: 

(3.76) Gxx < Gxy < Gyy < Hxx < ■ ■ ■ < Qyy < Gx < ■ ■ ■ . 
As a final general example including indices we have the inequalities 

m m 

(3-77) Hl^,^ < < GW,,,, <J2g' M,,, < Hi? H^^,^, 

k=l fc=l 
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extracted from (II) of Theorem 1.7 (3). 
In the sequel, we shall call 

• terms of order monomials like G, H, L M, G H M; 

• terms of order 1 monomials like Gx, Gx M, Lx 0; 

• terms of order 2 monomials like Gxy, Lyy, Mxx (our terms of order two will 
always be linear), 

according to the top order partial derivatives. 

3.78. A mean of checking intuitively the validity of partial differential relations. 

Before replacing (3.64) in the five remaining identities (3.65)2, (3.65)3, (3.65)4, 
(3.65)5 and (3.65)6, let us observe that if we assume that j ^ l\ in (3.69), then all 
the terms involving 6 vanish, so that we obtain the nontrivial partial differential equa- 
tions: 

(3.79) = -2G^^^+Hl^ + 2 5^ L^, ^1 

k k 

for j 7^ li. Here, we have underlined the first order terms plainly, in order to dis- 
tinguish them from the terms of order zero. These equations coincides with (I) of 
Theorem 1.7 (3), again specialized with j ^ h. Importantly, we notice that the choice 
of indices j 7^ li is possible only if m. ^ 2. Thus, we have derived a subpart of (I) as 
a necessary condition for the point equivalence to Y^-^ = 0, j = 1, . . . , m ^ 2. These 
first order equations show at once that there is a strong difference with the case m = 1. 

How can we confirm (at least informally) that the functions G^ , Hl^ and given 
by (3.35) in terms of X and do indeed satisfy these equations for j 7^/1? Dropping 
the zero order terms in (3.79) above, we obtain an approximated equation 

(3.80) Q^-2Gl,^+Hl^. 

Here, the sign = precisely means: "modulo zero order terms". We claim that this 
approximated equation is a consequence of the existence of X, Y'j . 

Indeed, according to the approximation (3.58), together with the definition (3.35) of 
the functions G^ and Hf^ , we have 



(3.81) 



Q3 = ^ _yj 



XX 



H = -2 , = -2 Y' . . 

1 1 1 xy 1 

Differentiation of the first line with respect to and of the second line with respect to 
X yields: 

(3.82) G^. ^-y^' and if/ ,^-2F^, , 

V / yii xxy'-i i-ii^ xy'-ix' 

SO that we indeed have = — 2 G^,^ + if/^ ^, approximatively and modulo the deriva- 
tives of order 0, 1 and 2 of the functions X, YK 

Similar verifications have been effected constantly in our manuscript in order to con- 
trol the truth of the formal computations that we shall expose until the end of Section 4. 

3.83. Continuation. From now on and up to the end of Section 4, the hardest compu- 
tational core of the proof may — at last — be developed. Further amazing computa- 
tional obstacles will be encountered. 
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Replacing plainly (3.64) in (3.65)2, we get: 



k k 



(3.84) 



+ 



+eB^^ + ^^'.0°)- 

Developing the products and ordering each monomial, we get: 

^ ® ^ ® 

+ -5/ l;^ , - + -(5/ e^e'i + 1 y //f Lf , - 

(3.85) 4 '^''i m 4 4 '2 Q 2 ^ '^''^ 

I — I k |Y| 

^ ^ ^ ^ 

--Lf, 6° +^5/L;^e° +^5/L;%e° +^5/e°e'i + 
^ ® 



We simplify according to our general principles and we reorganize the equality be- 
tween the first two lines of (3.84) and (3.85) so as to put all terms 6^; in the left-hand 
side of the equality and to put all remaining terms in the right-hand side, respecting the 
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(3.86) 



\ Si el^ + 1 si el^ - 1 si e?, 



1 



-- 

2 hd^ 



2 "'i ^y'^ 



k k 
k k 

+ \SIQ'Q^K 



Here, we have underlined plainly the four first order terms appearing in the second 
line. 

Next, replacing plainly (3.64) in (3.65)3, we get: 



(3.87) 



n 



h,l2,y''i 2 '2,/2,y'3 2 '2 h,h,y'^ 2 y^ 2 '2 y'a 

+ ^/l,/3,y'2 2 '1 '3,^3^2 2 '3 'l.'lV2 2 '1 S/'2 2 !^'2 " 



-M, 



2 ^3 2 '3 



E 



-Li,k + 2 



2 '3,^3 2 '3 



+ Ml,,, -[-^K + l ^'2 ®°) + E (-^' ,^3 + 1 € 
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Developing the products and ordering each monomial, we get: 
(3.88) 

^ ^ Jl h 



13 



■© 



1 



1 



1 



-® 



1 



4 "(2 "^i3,/3Q 4 Z3 



-® 



4 ^2 hi'l 



-® 



2 '3,h 



-® 



4 ^3 h,h 



10 



4 '1 '3, '3 . 



4 ^3 



-® 



4 '1 



2 '3,^2 



4 ^3 ^2,^2 



12 



4 '2 ^3,^3 



4 '3 



18 



4 '2 



-® 



1 2 I 1 15 I. I — I u I — I 



9 ^/l,/3 9 ^^2 5Z ^t«3 ~ ~ 



^«1,«3 



14 



i± — [3j ^± (d) i± Cs) ^± 



© ^ CD ^ ® ^ — 19| 



+ 1 K 



G) 



^ fi-i ^ -(k) ^ 



17 



4 



0/3 0/2 



11 



- si l!^ , e'l + - si e'^ e'^ + - si e'^ e'^ . 



We simplify and we reorganize the equality between the second and third lines of (3.88) 
and (3.85) so as to put all terms 0^ in the left-hand side of the equality and to put all 
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- si e'i, - - 5/ e\ + - si e\ - - si q\ 

o ll y'.i o (1 j/'2 o h y'3 n 13 y' 



(3.89) 



2 

— -I- i A-? —Ia-? I 

+ L\\ I. - -si L\\ . + 



Next, replacing plainly (3.64) in (3.65)4, we get: 



(nS,o)., - {Km) 



X 



2 ^''^ + 2 



E ^4,. - z E - 1 E 0' + 1 ^S,.. 0° + 1 0° 0"- 



Reorganizing the equality so as to put the terms 9^; and 6^, alone in the left-hand side, 
we get: 

(3.91) + E ^'i.'^ - i E - ^ E 0'+ 
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Next, replacing plainly (3.64) in (3.65)5, we get: 



^ I ^ I 

k k 

k ^ 



I — I I — I 1, I — I 



k VT\ k m ^ H] 



^ 

+ ^ l!i e'2 + ^ + ^ e'^ e'^ - - V i/,'; m,, ^ + 
^^i:^^ © © i © 2^ ^ " ^ 

^ — — — m 



Multiplying by —2 and reorganizing the equality, we get: 



-L 



h 



(3.93) 



k k 



CHARACTERIZATION OF NEWTONIAN FREE PARTICLES 
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(3.94) 



k k 
= -Mluh {\ ^iL^a + \ 0'^) - E ^^'3,. [-Lti2 + \ € Ltl2 + 

Developing the products and ordering each monomial, we get: 

= -l^th Mh,h- \ Mi,,i2 e'^ + J2 Li,i2 Mi,,k - I L\l,^ M,3,,, 



-® 



© 



(3.95) 



A (d) A © A © A © 

® _k ^ © © 

2 ® 2 ® 

Simplifying, we obtain the family (IV) in the statement of Theorem 1.7 (3): 



(3.96) 



= M 



M, 



3.97. Solving 6°, 6°,^, 6^^ and Q''^i,. It is now easy to solve all first order partial 

derivatives of the functions 0° and 0'. Equation (3.93) already provides the solution 
for . We state the result as an independent proposition. 

Proposition 3.98. As a consequence of the six families of equations (3.69), (3.86), 
(3.89), (3.91), (3.93) and (3.96) the first order derivatives 0°, 0°,^, 0^^ and Q^^^, of 
the principal unknowns are given by: 



(3.99) 



0*- 



2 E ^t>^ - E - 1 E Ht- 

k k k 
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(3.100) 



f)0 - 1 rh _ ± ttI 



2 ti,ti,a) 2 



fc k 

+ 3 2^ ^'i ^'i.*^ ~ 2 ^ ^h^k,k - 2 2^ + 
fc fc 

+ lLi^ 0O + l@o@h_ 



(3.101) 



(3.102) 



2 1 



2 '^''i 2 



9'i 



h 



k k 



We notice that the right-hand side of (3.99) should be independent of h; this phe- 
nomenon will be explained in a while. 

Proof. For 6° in (3.99), it suffices to put j := h in (3.69). 

To obtain 6°;^, we put j := I2 and I2 := h in (3.86), which yields: 



X 2 y'l 



1 W'l -L 
'2 



(3.103) 



fc 

fc A: 

4 '1 4 'i.'i 4 

We may easily solve 6°,^ and Q''^ thanks to this equation (3.103) and thanks to (3.91): 
indeed, to obtain (3.100), it suffices to compute |- (3.91) + |- (3.103); to obtain (3.101), 
it suffices to compute | ■ (3.91) + | ■ (3.103). Finally, (3.102) is a copy of (3.93). This 



completes the proof. 



□ 
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3.104. Appearance of the crucial four families of first order partial differential 
relations (I), (II), (III) and (IV) of Theorem 1.7 (3). However, in solving 9°, 9°,^, 

9^1 and Q^^^^ from our six families of equations (3.69), (3.86), (3.89), (3.91), (3.93) 
and (3.96), only a subpart of these equations has been used. We notice that the two 
families of equations (3.91) and (3.93) have been used completely and that the family 
of equations (3.96), which does not involve 9, coincides precisely with the system (IV) 
ofTheoreml.7 (3). To insure that 9°, 9°;^, 9^^ andQ^^i^ as written in Proposition 3.98 
are true solutions, it is necessary and sufficient that they satisfy the remaining equa- 
tions. Thus, we have to replace these solutions (3.99), (3.100), (3.101) and (3.102) in 
the three remaining families (3.69), (3.86) and (3.89). 

Firstly, let us insert inside (3.69) the value of 9°, given by the equation (3.99), in 
which the index li is replaced in advance by an arbitrary index I2. We get: 



+ 2 E 4,, L\lk -^IY.G' Ll 



-© 



(3-105) + ^ _ 1 ^ + 1 ^ Hi Hi^- 



5^ y G^Q^ + (5/ y G''Q^ + - 5^ 9° 9° - - 5/ 9° 9° . 
® ® ® ® 



We simplify, which yields the family (I) of partial differential relations of Theo- 
rem 1.7 (3): 



(3.106) 



k k 



Secondly, let us insert inside (3.86) the values of 9^^, 9^^ given by (3.101) and 
the value of 9°;^ given by (3.100). We place all the terms in the right-hand side of 
the equality and we place the first order terms in the beginning (first three lines just 
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below). We obtain: 
(3.107) 

= --H^ +Lj 
+ - 6^ H^' 



h,l2,x\ 



61 L 



h 



[T| 2 'i '2.^2,x u| 2 



h h,li,xi 







+ 



- 5/ l!^ 



Q <2 Zi.ti'ir 



Li' 



+ 



1 



^_h_^2hr^^ Q h l2,y^2 







+ 



k 

k,k 



12 



13 



-© 



i 5i y Hi +\ 5/ lY , e° + ^ 5/ 6° G'l - 

4 t2 'l 4 '2 tl.'l ^ 4 '2 ^ 

^ (b) ® ® 



l2,k 



10 



14 



+ 



15 



^ © ^ © ® ® 

-l^iG'^M,,, -l6iT.G'M,,,, +lsiY.Hi'Ll,^ -U^T.^^,!^ + 
I'^l k I — \ k I — i fc I — I 

I 1 K 1^ K K 

^ © ^ © ® ® 

' — ' ^ [To] - [T4] *i TTF 

-^^E^4^ -i^^E^e' +^^^t20Vi'^®°®'^o- 
t © 1 © ^ ® 
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Simplifying and ordering, we obtain the family (II) of partial differential relations of 

Theorem 1.7 (3): 

(3.108) 

- W + - I 1 d rrh I 

u - 2 ^i-y'^ 6 3 'iV^ 

1 . 9 . 



k 

k k k 

k k 



2 -(2 / u K 2 '12 / J '1 ~^h,k' 

k k 



Thirdly, let us insert inside (3.89) the values of e'^i,, of Ojjl^, of Q^^^^^ and of Oj^^^ 
given by (3.102). We place all the terms in the right-hand side of the equality and we 
place the first order terms in the beginning (first four lines just below). We obtain: 



0=-S^ , , - Ml, (5/ , , - 5/ Ml, I, 

2 ^1 ^2,/2,y'3 Q h i-2,i.3,x ^ 2 li,li,y^3 Q h 'i.'s.a^ jyj 



2*1/3^34/2 h l-3,l2,x ^ 3j3Zi^l4/2 ;3 ti,(2,3: rji 

(3.109) ^ 7^ 

4- r-' — r-' -I- - P — -/i-? r'2 1 



1 



1 1 ^/3 _ 1^^- H2, 



i 



+ 



11 



^ rk T k ^ xi T k T k 

+ 9 % 2^ -'"«l:^2 ^k,k - n"l2 2-^ ^h,h ^k,k + 



-® 



-© 



+ hiL\], @^^ -l5id],e'- +i5/L!^e'^ -lsldi,e^^ + 



4 '2 



4 '3 



^ 4 '2 ^3,^3 ^ 4 '3 ^2,/2 



-© 



2 ^'^ © 2 '3 'l^'2 ^ ^^2 Q ^J3 Q 
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^ ® ^ © 

© © ^ k ^ 

© 

2 " ^ © © 

® — © 



10 



^ © ^ © ^ k 



^ ® ^ © 

+ 1 ^ti. 0" + 1 ^iL^. 0" - 1 5: 4,.. 0' + 

© © ^ k Q 

^ © ^ © 

^ ^ ^ © 

+ i ^t. 0' + 1 4,. 0' - E 4... 0^= + 

^ © © ^ k ^ 

© 

+ - 5/ i, e° +\ 6i . 
2 " ^ © ^^f © 

Simplifying and ordering, we obtain the family (III) of partial differential relations of 
Theorem 1.7 (3): 



1 , . 1 
2 



(3.110) 



+ ^ '^/x E ^'3 ^^..^^ -l^lH M,,,,+ 

k k 

k k 

+ E ^h,h ^h,k ^ E 4,/2 ^h,k- 
k k 

3.111. Arguments for the proof of Theorem 1.7 (3): necessity and sufficiency of 
(I), (II), (III), (IV). Let us summarize the implications that have been established so 
far, from the beginning of Section 3. Recall that m ^ 2. 
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• There exist functions X, of (x, y) transforming the system 
vix = F^{x,y,yx), j = l,...,m, to the free particle system Y^-^ = 0, 
j = 1, . . . ,m. 

• There exist functions ^j-^i^ of ix,y), ^ i,h,h ^ satisfying the first 
auxiliary system (3.38) of partial differential equations. 

• There exist (principal unknowns) functions 6°, 0^ satisfying the six families 
of partial differential equations (3.69), (3.86), (3.89), (3.91), (3.93) and (3.96). 

• The functions , H^^, Lj^ and satisfy the four families of partial dif- 
ferential equations (I), (II), (III) and (IV) of Theorem 1.7 (3). 



The four families of first order partial differential equations (3.99), (3.100), (3.101) 
and (3.102) satisfied by the principal unknowns will be called the second auxiliary 
system. It is a complete system. 

To achieve the proof of Theorem 1.7 (3), we have to establish the reverse implica- 
tions. More precisely: 



• Some given functions G^, Hj^, Lj^ ^^ = L\^ i^ and Mi^^i^ = Mi^^i^ of {x,y) 
satisfy the four families of partial differential equations (I), (II), (III) and (IV) 
of Theorem 1.7 (3), or equivalently, the partial differential equations (3.106), 
(3.108), (3.110) and (3.96). 

• There exist functions 6°, satisfying the second auxiliary system (3.99), 
(3.100), (3.101) and (3.102). 

• These solution functions 9°, satisfy the six families of partial differential 
equations (3.69), (3.86), (3.89), (3.91), (3.93) and (3.96). 

• There exist functions ^j^i^ of {x,y), ^ j,h,h ^ satisfying the first 
auxiliary system (3.38) of partial differential equations. 

• There exist functions X, Y'j of (x, y) transforming the system 
yl^ = F^{x,y,yx), j = 1, . . . , m, to the free particle system Y^;,^ = 0, 
j = 1, . . . ,m. 



The above last three implications have been already implicitely established in the 
preceding paragraphs, as may be checked by inspecting Lemma 3.40 and the formal 
computations after §3.62. 

Thus, it remains only to establish the first implication in the above reverse list. Since 
the second auxiliary system (3.99), (3.100), (3.101) and (3.102) is complete and of 
first order, a necessary and sufficient condition for the existence of solutions follows 
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by writing out the following four families of cross-differentiations: 



(3.112) 



9<i 



yh 



9'i 



V y 



In the hardest techical part of this paper (Section 4 below), we verify that these four 
families of compatibility conditions are a consequence of (I), (II), (III) and (IV). For 
reasons of space, we shall in fact only study the first family of compatibility condi- 
tions, i.e. the first line of (3.112). In the our manuscript, we have treated the remaining 
three families of compatibility conditions similarly and completely, up to the very end 
of every branch of the coral tree of computations. However, we would like to men- 
tion that typesetting the remaining three cases would add at least fifty pages of Latex 
to Section 4. Thus, we prefer to expose thoroughly the treatment of the first family 
of compatibility conditions, explaining implicitely how to guess the treatment of the 
remaining three. 



§4. Compatibility conditions for the second auxiliary system 

So, we have to develope the first line of (3.112): we replace 0° by its expres- 
sion (3.99), we differentiate it with respect to ?/'\ we replace 0°,^ by its expres- 
sion (3.100), we differentiate it with respect to x and we substract. We get: 
(4.1) 

+ 2 E Gl.. Liu + 2Y^G^ - E L\u - E i^^,,,, - 

k k k k 

k k 

+ 9090, ■ 



y 

2 l\,l\,xx 2 /i.J/'ix 



3 ^ X .1,^ 3 

k k 



+ 3 ^k,x ^h,h ^ ^ ~ 3 ^ 'l'^ 'i''^ ~ 3 2^ '1 li,k,x+ 

fc f\> Jij f\> 

+ 9 E ^k,k + 9 E Lk,k,x + 9 E ^h,x + 9 E 



2 / _j Ll,X — K,K ' 2 / ^ if K^K^X • 2 / ^ l-l',-^ 2 

Aj Aj 

^ W ©° - ^ ^w. - ^ eS 9'^ - i 90 9^. 
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Here, we underline twice the second order terms. Also, we have underlined once the 
six terms: ©^ii, ^yh^ ©x' ©2' ©S ©x- They must be replaced by their values 

given in (3.99), (3.100), (3.101) and (3.102). In this replacement, some double sums 
appear. As before, we use the first index A; = 1, . . . , m for single summation and then 
the second index p = 1, . . . ,m for double summation. Finally, we put all the second 
order terms in the first line, not disturbing the order of appearance of the 73 remaining 
terms. We get: 
(4.2) 

I ^ I 2 ; 



^ - ^^ ^ -A 

^ \^ n'' fr'i ^ n'' w'l \" r^fc cvfc , 



^yh ^k,k - l^'^ ^k,k,y'i ~ 
^ © 



10 



11 



18 



fc p 



19 



1 \ ^ jk jh I \ " \ " ^fc rP rp 1 \ " f-ik rk cdi 

2 ^k,k^l,,l, +2^2^^ ^k,h^P,P ~ 2 2^ ^k,k^ ~ 
^ © Jl E ^ ^ © 



jE^'^iLh©' +EEg'^L0' -j:G^M,,,eo -i J^G'^e'^e'^ + 



-© 



A; p 



-© 



-© 



-® 



3 '"'^'^ ® 3_Ji:]^® 



1 



-© 3 k 



-© 



-® 



^E^u^t0° -^E^u0'0° +^e°e°e'^ 
(D © ® ' 



A A 



17 



-0 



14 



13 



+ ^ E ^h,^ ^k,k + n E ^k,k,x + ^ E ^h,x 



-0 



15 



-© 
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q Yl ^k,y^ + fi X] ^k,k,x + o X] ^k,k + 



12 



+ -\EE^tH',Ll, 



k p 



21 



k p 



k p 



25 



lEE^i.H'^^'p. -zEE^^^r©^ +iE^u4.0° + 



k p 



26 



fc p 



-CD 



E^fc T-ii rZi I 1 \ ^ ^fc fh jk , -'- \ ^ rrk Tjh jh i 

<^ ^h,h ^li,k +9 2^*^ + 7 2^ ^'i ^fc + 



22 



27 



ij^ ® ^ ® 

-El (7!\ 2 'i'^ r:^ 



EG^^I^.fe©'^ +jE^'^t©'' +iE^';^fc©'' +i^G'=e'^e'^ 



-® 



-© 



-® 



i © 

^ ^ © © ^ © 

^ ® ^ © ^ ^ 

As usual, all the terms underlined with the 15 roman alphabetic letters a, b, . . . , n, o 
appended vanish evidently. Furthermore, we claim that the eight terms underlined 
with the 8 Greek alphabetic letters a, (3, 7, 5, e, (, rj and 9 also vanish: 
(4.3) 

=?= - y G^' + + - y i^/; ^. 9^= - - hI^ ^ e'^+ 

/ V J/'i y'l 2 / ^ 1' 2 ' 

+ E E ^1.. ©' - E ^iL. ©'^ - i E E K Hie^ + \Y: K ^i^ ©'^ • 

k P k k P k 

Indeed, it suffices to observe that this identity coincides with 



(4.4) 



= ^E ©'((3-106)1,,=,; 



h-=h; l2-=h 
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Simplifying then (4.2), we get the explicit formulation of the first family of compati- 
bility conditions for the second auxiliary system: 

4 2 
n — ?— 9 r^'i 4- _ r'l — 



3 ^'^ 3 

fc 



(4.5) 



k k 
1 trii r«i 1 \^ 7J-A.' r/i , ^ fjk rk , 

k k k 

~ 3 ^k,yk^h- 

k 

3 ^h,k,x^h + 3 2^ ^h,h,x^k +32^ ^k,k,xJ^h + 

k k k 



- E E ^'.p + \Y.g' K Mk,k + ^ E E 

k p k k p 

~ E ^'L^i + E E '^'^ -^Mi 
k k P 

4 5: E ^1. + i E E K 

k P k P 



k P k 

We can now state the main technical lemma of this section and of this paper. 

Lemma 4.6. The second order partial differential relations (4.5) hold true for I = 
1, . . . , m, and they are a consequence, by differentiations and by linear combinations, 
of the fundamental first order partial differential equations (3.106), (3.108), (3.110) 
and (3.96). 

4.7. Reconstitution of the appropriate linear combinations. The remaining of Sec- 
tion 4 is entirely devoted to the proof of this statement. From the manual computational 
point of view, the difficulty of the task is due to the fact that one has to manipulate 
formal expressions having from 10 to 50 terms. So the real question is: how can we re- 
constitute the linear combinations and the differentiations which lead to the goal (4.5) 
from the data O-'^Oe), (3.108), (3.110) and 0.96)1. 

The main trick is to first neglect the first order and the zero order terms in the 
goal (4.5). Using the symbol "=" to denote "modulo first order and the zero order 
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terms", we formulate the following sub-goal: 

3 li,xyh 2 



(4 8) =?= -2 G'^. . + - i/^ - - ^1 



for /i = 1, . . . , m. Before estabilishing that these partial differential relations are a con- 
sequence of the data (3.106), (3.108), (3.110) and (3.96) (written with a similar sign 
=), let us check that they are a consequence of the existence of the change of coordi- 
nates (x, y) (— (X, Y) (however, recall that, in establishing the reverse implications 
of §3.111, we still do not know that such a change of coordinates really exists); this 
will confirm the coherence and the validity of our computations. Importantly, we have 
been able to achieve systematic corrections of our computations by always checking 
them alongside with the existence of the change of coordinates (x, y) i— > (X, Y). 
Coming back to the definition (3.35) and to the approximation (3.58), we have: 

(^li _ rih ^ y^'i 

^ — XX XXI 

(4.9) K = -2nlVi + °- - -2 ^xV. + 

Differentiating the first two lines with respect to y'^ and the third line with respect to x, 
and replacing the sign = by the sign = (in a non-rigorous way, this corresponds essen- 
tially to neglecting the derivatives of order 0, 1, 2 and 3 of X, Y^ and to neglecting the 
difference between the Jacobian matrix of the transformation and the identity matrix), 
we get: 

qIi _ _Y^i 

yHyH xxy^iy^i ' 



(4 10) h=-'^Y^\, h+^xxxvh 

\T.±\jj li,xy'-i xxxy ^ 



l i,li,xx — y'lyhxx xy'-ixx- 

3^ 



result: 

4 2 

U ^^yhyh 3 ^h,xyh 3 ^h,h,x 



Hence the linear combination -2 • (4.10)i + | ■ (4.10)2 - 1(4.10)3 yields the desired 



(4.11) =9y'i +lv , + ^ y'l --X , 

- ^xxyhyh ^ 3 xyhxyh ^~^ 3 ^xxxy'l ^T 3 ^hyh^^ ^ 3 ^xy'lxx ^ 

= 0, indeed! 

Thanks to this straightforward computation, we guess that the approximate partial 
differential relations (4.8) are a consequence of the approximate relations (3.106), 
(3.108), (3.110) and (3.96), namely: 



(3. 106)-°^^ : 0^-2 G^,^ + 2 5/^ G[% +Hl^- 6^ H\1 



2 

yii ' ^'^h^yh ' '^'^h,x "h '^'^l2,x^ 
1 1 1 



(3.108)"^°^ : = — H^, u+- 5/ H]^ + - 5/ H]^ u 

1 
3 



(4-12) , ri _lxir«2 Th 

^h,l2,x o"h^l2,l2,x n"h-^h,li,x^ 



(3.110)-°^^ : = L[^^^^^,3 - L^^^^_^, + 5^ M,,,,, - 5^ M,,,3„ 
(3.96)-d : = M,^,,^,,.3 - M,^,,3,,., 
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Here, the sign = means "modulo zero order terms'". Before preceding further, recall 
the correspondence between partial differential relations: 

(I) = (3.106), 

(II) = (3.108), 

*• ■ ^ (III) = (3.110), 

(IV) = (3.96). 

However, these couples of equivalent identities are written slightly differently, as may 
be read by comparison. To fix ideas and to facilitate the eye-checking of our subsequent 
computations, we shall only use and refer to the exact writing o/(3.106), o/(3.108), 
o/(3.110) and of 0.96). 

4.14. Construction of a guide. So we want to show that the approximate relation (4.8) 
is a consequence, by differentiations and by linear combinations, of the approximate 
identities (4.12). The interest of working with approximate identities is that formal 
computations are lightened substantially. After having discovered which linear com- 
binations and which differentiations are appropriate, i.e. after having constructed a 
"guide", in §4.22 below, we shall write down the complete computations, including all 
zero order terms, following our guide. 

We shall use two indices li and I2 with 1 ^ /i, /2 ^ m and, crucially, h 7^ h- Again, 
the assumption m ^ 2 is used strongly. 

Firstly, put j := li in (3.106)™°'^ with I2 7^ h and differentiate with respect to ?/'^: 

(4.15) = -2G'i. . +2G\ . u-H]^ I,. 
^ ' y^y^ y^y^ h,xy'i t2,xj;'i 

Secondly, put j := I2 in (3.106)™°*^ with li 7^ I2 and differentiate with respect to y'''^: 

(4.16) = -2G\ i^+H]^ u- 
^ ' y^y^ h,xy'2 

Thirdly, put j := I2 in (3.108)™°^ with li 7^ I2 and differentiate with respect to x: 

Fourthly, put j := h in (3.108)™°'^ with I2 ^ h and differentiate with respect to x: 

^^■^^^ ^ = ~2 ^h,y'2x + 6 ^h,y'2x + ^h,l2,xx " 3 ^h,l2,xx- 

Fithly, permute the indices (/i, I2) ^ (hJi)' 

^^-^'^^ ^ ^ ~2 ^h,y'^x + Q ^h,y'ix + ^h,h,xx - 3 ^h,h,xx- 

Finally, compute the linear combination (4.15) + (4.16) + 2 ■ (4.17) - 2 ■ (4.19): 

n — 9 Wi _|_ 9 /^'2 _|_ zrh zrh 



2G\ . +Hj^ 

y^y^r^ h 



(4.20) , 2 , , 4 , 

— H -I- — H -4- 2 T — — T -\- 

ZIlMjllQ 3 h,luxyh ^ ^ ^^h,l2,xx ^ 3 ^h,h,xx ^-^ 

1 2 

l2,xy'i ^^ 3 li,xyh ^ l2,li,xx (^'^ 3 luluxx ^ 
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We indeed get the desired approximate identity: 

(4.21) ^^-'^G'yHyH+\K,^yH-\L\\,^,..- 

4.22. Complete computation. Now that the guide is constructed, we can achieve the 
complete computations. 

Firstly, put j := li in (3.106) with I2 7^ h and differentiate with respect to y'^: 
(4.23) 

0^-2 G\ , + 2 G\ , + Hj' , - Hj' , + 

y^y^ y^y^ h,^y^ h,xy^ 



Jij A/ 

--Y H!r . f!^ - - V //f i^fs + - y //f . 4^ + - y ^f; , . 



2 ^ '1:3/^ 2 ■'^ ^ ^''f 2 2 

Ay A/ Aj ^ 

Secondly, put j := /2 in (3.106) with li ^ I2 and differentiate with respect to y'^: 
(4.24) 

y'lyh h,xyh 
A A A A 

Thirdly, put j := I2 in (3.108) with /i 7^ /2 and differentiate with respect to x: 
(4.25) 

11 2 

2 li,y'2x^ ^^h,yhx^ ^h,h,xx ^^h,h,xx^ 

+ &i Mi.^i, + G'i Mi,j,,x - ^ G'^ Mi,j, - ^ Mi,j,,x- 



fc fc fe k 

I 1 trfc ri2 I fj-fc rZ2 i ^ fj-'l T _L ^ IT^l 

k k ' 

^ \^ ^ \^ r'l 



A^ A A A? 



Fourthly, put j := /i in (3.108) with /2 7^ /i, differentiate with respect to x and permute 

the indices [h-ih) ^ (hji)'- 

(4.26) 

^ ^ 2 l2,y^ix g h,yhx ^l2,ll,xx 3 

+ M,,,z, + G'2 Mi,,;,,, - ^ G^i Mz,,i, - ^ G'l Mi,j,^x+ 



+ IEg'x Ml,,, + Iyg' Mi,,,,x - 1 Y <x -lE^k lU,x+ 

J'C A A 

1 1 Yji^ 1 TjW rk ^ ttW r A; 

+ 2 2^ '2,2: -^ii,fc + 2 2^ '2 ^h,k,x + e 2^ + g 2^ -'^fc ^ii,iux- 

k f\> k I\j 

^ \^ r'l ^ w'^' r'l 
~ R 2^ ^'i,^' ~ fi 2^ -"'i ^U, 



g / ^ ll,X~li,K g / ^ tl "tTfc,X" 

k k 
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Finally, compute the linear combination (4.23) + (4.24) + 2 ■ (4.25) - 2 ■ (4.26): 
(4.27) 

I ^ I 2 / 

= - '^'^yhyh + 3 ^luxyh - 3 ^l\,h,xx - 

k 

+ E -^'t^: -^la.fc 3 E ^k,x ^h,li ~ 3 E ^h,x ^h,k ~ E ^l2,x ^h,k + 

l\y /C A> 

A; fc 
+ E ^\l,k,x + 3 E ^h,h,x ^fc ~ 3 E ^h,k,x ~ E -^'Lfc.i: -^'2 + 

J\/ 

+ 2 y , + 2 y _ 2 y l'^ , G*^- 

/ ^ h,k,y'^ / ^ h,k,y''2 / j l2,k,y'-i 



3 . 



In this partial differential relation, importantly, the second order terms are exactly the 
same as in our goal (4.5). Unfortunately, the first order and the zero order terms are 
not the same. 



4.28. Formulation of a new goal. Thus, in order to get rid of the second order ex- 
pression 2 G'^,^^,^ + I Hll^y,^ - I L\li^^^^, we substract: (4.5) - (4.27). In the result, 
we write the first order terms in a certain way, adapted in advance to our subsequent 
computations. For this substraction yielding (4.29) just below, we have not underlined 
the terms in (4.5) and in (4.27). However, they may be underlined with a pencil to 



56 



JOEL MERKER 



check that the result (4.29) is correct. We get: 

(4.29) 

=?= 

1 1 

\ ^ r<k jk I /-ih fh I \ ^ ifk rk ^ ttIi jh i 

- S'l + J/'i h,h + 2 ^ 'i'^ ''''' ~ 2 'i'^ 'i.^i"*" 
_fc k 

k k 

k k 

fxi Jc 

^ \ ^ J-r'2 fjfc \ ^ r(2 uk . 

~ 2 ^ '^Vi '2 + 2^ ^h,k,x^h + 

k k 

+ -Y Hf ,^Hi^--Y uHl^ + 
fc fc 



- E E Hi Ml,,, + ^ E ^^^.'^ + ^ E E < ^^.p- 

k p k k p 

~ E '^'^ ^iLfc + E E -^p.p ~ 3 E E ^k,k+ 

k k p k p 

+ ^ E E K Hi ^Iv - i E E Hi L% + 1 Ri H^^ L\l . 

k P k P k 

We have underlined plainly the first order terms appearing in lines 1, 2, 3, 4, 5, 6 and 
7. 

4.30. Reconstitution of the subgoal (4.29) from (3.106), from (3.108) and 
from (3.110). Now, it suffices to establish that the first order partial differential 
relations (4.29) fox 1 ^ li,l2 ^ m and I2 7^ h (crucial assumption) are a consequence 
of (3.106), of (3.108) and of (3.110) by linear combinations. The auxiliary index I2, 
which is absent in the goal (4.5), will disappear at the end. Differentiations will not be 
applied anymore. Also, the partial differential relations (3.96), which were not used 
above, will neither be used in the sequel. However, they are strongly used in the treat- 
ment of the remaining three compatibility conditions (3.112)2, (3.112)3 and (3.112)4, 
the detail of which we do not copy in the typesetted paper. Finally, the construction of 
a guide for the subgoal (4.29) may be guessed similarly as in §4.14 above. We shall 
provide the final computations directly, without any guide: they consists of the seven 
partial differential relations (4.32), (4.33), (4.35), (4.37), (4.39), (4.43) and (4.45) 
below. At the end, we shall make the addition (4.47) below, producing the desired 
subgoal (4.29) := (4.32) + (4.33) + (4.35) + (4.37) + (4.39) + (4.43) + (4.45), with the 
numerotation of terms corresponding to the order of appearance of the terms of (4.29), 
as usual. 
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Firstly, put j := k, li := li and I2 ■= h in (3.106): 
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(4.31) 



p p p p 



Apply the operator ^ J2k ki') the preceding equality, namely compute 
I 4. ■ (4.31). This yields: 

n \ ^ r^k jk , /-ih I ^ \ ^ ijk jk ^ Tjh jh i 

^ - ^y'l ^k,k + '^yh + 2 ^ 'i'^ ''''' ~ 2 'i'^' 

_fc fc 

(4.32) + E E - E ^Iw - i E E ^r, ^^p' 

k P P k P 



4 

p 



Secondly, apply the operator — ^(■) to (4.32), namely compute — Yl,k k 

(4.32). This yields: 

(4.33) 

/ J j/'i i2:fc j/'i '2i<i / > 'ii^; t2,fc ti,a; t2,ti 



2 E E 4. + 2 E lIp+\T.Y. K h'p 

k p 

2 E Hp ^h,li- 



2 

k p p k P 



Thirdly, put j := k, li := I2 and I2 := /i with /2 7^ /i in (3.106): 

= - 2G;,, + 2^f^ + Ht^^^ - 5l H\l^ + 

^"■'"^ + 2 E «^ ^t. - 2 E <p - \ E K + 1 K ■ 

p p p p 

Next, apply the operator L|J ;j(-) to (4.34), namely compute L|J ^ ■ (4.34). This 

yields: 

(4.35) 

= -2 y + 2 g\ l\\ + y ^ - hi^ + 



+ 2 E E «^ ^tp <k - 2 E ^S. ^i? ,.2 - ^ E E K ei 

fe p 

+ 9 E ^fi-^p ^«?,«2- 



2 

k P p k P 
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Fourthly, put j := h, h := k and I2 ■= h in (3.108): 
(4.36) ^ ^ ^ 

1 ^2 ^ T k 

p p p p 



Next, apply the operator — Hi^{-) to (4.36), namely compute — Hf^ ■ (4.36). 
This yields: 

(4.37) ^ ^ ^ 

k k 

Erh TTk _1_ ^ r'2 [7^2 I ^ \" rfc rjfc 

^k,l2,X'"h+ ^^l2,l2,X^h +3 2-^ ^k,k,X^h- 
k k 

- G'^ Hi Mk,i, + ^ G'2 M;,,^, + ^ E Mk,k- 

k k 

-l^G' H\l M,„, + ^ E E G'^' ^^.p + ^ E E Hi 4^ L^.,- 

p k p k p 

2 / / ^ '1 fe '2,P g / > ti P '2, '2 g / h h h,P 
k p p p 

-lEE Hl^U + ^ E E Hi^ir 

k P k P 



Fithly, put j := I2, h := k and I2 := h in (3.108): 



''^ TJ^'^ J- ^ a'2 I/'i _|_ r '2 ^ r'2 r'l _|_ 



(4.38) 



+ G'2 M,,,, - i 42 G'^ M,,,,, + ^ 4^ E MiuP - ^ E + 

p p 

+ IEhI ^Ip + E ^p^ ^l.. - e k a 
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Next, apply the operator J2k ^hi') (4-38), namely compute J2k ' (4-38). This 
yields: 
(4.39) 



= -- V , Hi + - I + y Ht 



3 ^h,h,X^l2^ 



k p 

- ^ E E '^f. ff? ^Im + i E E < ''J '^i;,, ^\T."'i ■ 

A: p 

- - y „. 



k p 



Sixthly, we form the expression: 



(4.40) 



(3.108) 



\j:=k; h:=h; l2-=l2 



- (3.108)1 



j:=k; h--=h\ h-=h 



Writing term by term the substractions, we get: 
(4.41) 







1 



+ - Hi u + - <5f, h]^ 
o h,y'-^\ — I ft '1 I' 



1 



St Hi' 



' ~ H, i„ -| jn, I, -| 111, £1, /„ — — n, 







+ 



+ l5tH'- 



1 



k Tjh 



- St H 



3 3 ^2^2 







1 



1 



^ rfc r'2 _i_ ^ A*^ 







L 



l2,h,x 



k jh 







+ 



10 



I 6l Y: M,,, + 1 M,,,, - \ Y: Hi, Ll,^ +\YHi LI,. + 

[E © © 



12 



13 



14 



^/';E^r2^L +^tE^f.^L 



15 



16 



14 



^^/';E^^^L2 -Hek^Ip +\^tEKLtp 



13 



16 



15 
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Simplifying, we get: 
(4.42) 



p p 

+ \ E ^f. ^Iv <v - \ € E + 1 € E 

p p p p 

+ l<EKLtp-l€EKLL- 



Next, apply the operator - Y,k H^i-) to (4.42), namely compute - Y.k ' (4.42). 
This yields: 
(4.43) 



2 ^-^ h,y'^ «• 2 '2,s/'i K g /2,j/'2 /i 

k k 

_ uh Tjh _^ T^2 Tjh I ''^ r'l _ 
g -"ii^i ^2 3 ^h,l2,x ^ 3 ^hA,x "-h 



p 

+ H\i - ^ E E K ^Ip + ^ E E ^f. 

P k P k P 

+ -y Hf' Hi' Lf , --Y Hi' H\}L^ , --Y Hi' Hf L\' + 

g / ^ (i P h,h Q / y h P h,h g / > ti t2 (2,p 

p p p 

+ - V if?^ 



Seventhly, put j := I2, h := k, I2 := /2 and /s := /i in (3.108): 



Mk,h,x+ 



(4.44) 



p p 
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and then apply the operator 2 J2k G'^'(-): 



(4.45) 



fc fc p 

- G'^ i/f^ M,,,, - E E Hi M,,,, + 2 5: L^,^ Lj^,^- 



k p 



k p 



k p 



Finally, achieve the addition 



(4.46) 



(4.32) + (4.33) + (4.35) + (4.37) + (4.39) + (4.43) + (4.45). 



We copy these seven formal expression, we underline the vanishing terms and we 
number the remaining terms so as to respect the order of appearance of the terms of 
the subgoal (4.29): 
(4.47) 

^ - S'l ^.'^ +%'irr|+ 2 ^ 'i'^ ''''' ~ 2 'i'^ 1 

-h . □ ^ ^ . ^ 

+ E E ^h,P ~ E ^h,h ^h,P ~ 7 E E Hp ^k,k + 



k P 



24 



23 



k P 



27 



28 



I o \ ^ r^k rh o J h \ ^ rrk rh I uh T h 



-0 



-2EE^'^t^Lp +2Eg'</.^L +^EE^f.^P^: 



k p 



\ E K ^ti. - 
© 



k p 



-CD 



2 y G^, . +2 G'\ l!% + \" if,'^ ^ l!^ - //'i , + 

+ 2 E E Lf,, Ll, -2YG^ lIp ^ti. - ^ E E K H'p Ltk + 
© ® -^-^ ® 



k P 
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I ^ \ ^ ijh ijk -'- Tjh Tjh ^ \ ^ Tjk Tjk 

k nr\ k 



10 



1 2 
_fc vjn ® Vfc 



12 



Y^G'^K^k^i, +\g'-h\IMi,^i, +lY,G^K^iKk - 



^ © ^ ^ ^ 



-© 



^EE^^^^^^P -^E^^^^L. +^E<^f2^t 
© 



k P 



-© 



-© 



k P 



25 



26 



6 " © . 



13 



E ^k,h,x 



14 



--Li\ H]^ + 

o h,h,x h „ 
^2 



-Q 



p 



k p 



-© 



k p 



i V H]' Hf lY + 

g / ^ t2 (i h,p 

© 



Ij?] lie] ~ 



© 



6 '2,»/'2 '1 



1 



i_r'i Tjl'2 ^ u'2 I ""^ r'l W 



6 



•© 3 



© 



h 



© 



-© 



© 



+ ^Eg^<^^'..p -\Y.Y.H'kKLlp +\Y.Y. H'k Ki^lv + 



I © © 

g / ^ ti P h,h g / ^ '2 P tli'l g / ^ tl '2 t2,P 

^ © © © 



k p 



-® 



+ ^E<^^t + 
© 
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17 



18 



19 



+ Y,G''HllMi,,i, -Y,G^H\lMk,i, +Y,G'^HlMi,,, - 
^ © ^ © ^ © 



-© 



k p 



20 



k P 



k p 



-© 



In conclusion, there is exact coincidence with the subgoal (4.29). The proof that the 
first family (3.112)i of compatibility conditions of the second auxiliary system (3.99), 
(3.100), (3.101) and (3.102) are a consequence of (I), (II), (III) and (IV) of Theo- 
rem 1 .7 (3) is complete. Granted that the treatment of the other three families of com- 
patibility conditions (3.112)2, (3.112)3 and (3.112)4 is similar (and as well painful), 
we consider that the proof of the equivalence between (1) and (3) in Theorem 1.7 is 
complete, now. □ 



§5. General form of the point transformation 
of the free particle system 

This section is devoted to the exposition of a complete proof of Lemma 3.32. To start 
with, we must develope the fundamental equations (3.10), for j = 1, . . . , m. Recalling 
that the total differentiation operator is given hy D = ^ + YlTi=i V^x ' + YmI=i V^xx ' 
-T^, we compute first 

(5:i) 



DDX =D 



m. 

Vx Vx ' -^y'lyh + E^ yxx ' -^y'l J 



h=l 



lj_=l l2 = l 



and 
(5.2) 



DDY^ =D 



yL 



m. m 

E y'^' ' ^xyh + E E y'^ ■ -^iiy^a + E y'^'^ ■ -^j/i 
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Now, we can develope the equation = —DY^ ■ DDX + DX ■ DDY^ , which yields 



h=i 

m m 



Xxx + 2 ^ y^^ ■ X^yii - 



(5.3) 



h=l l2 = l h = l 



+ 



Xx + Y,y'^"- ^y'- 



m m 



/l = l ^2 = 1 



Zl=l 



-XxxYi + YLXx+ 



^xx J- X ^ ^ XX 

m 



h=i 



m m 



2X,y.,Y^ + 2Y^^,^Xx- 

— X Y^ -\- Y^ X , 

^xx 1 yi^ ^ 1 XX ^s/'i 

— Xyl^yl2 Y^ + Yyl^yl^ 



-2 X„..i2 Y\^ + 2 Y'' ,„ X.i-i^ 



xy'^ -yi 



m m m 

il = l l2 = l l3 = l 
m 

m m 

+ J2T.y'^'^y'^'-[-^y'^K^^ + K^^^y 

ll=l «2 = 1 



+ 



The goal is to show that after solving these m equations for j = 1, . . . , m with re- 
spect to the yl^, I = 1, . . . , m, one obtains the expression (3.33) of Lemma 3.32, or 



CHARACTERIZATION OF NEWTONIAN FREE PARTICLES 



65 



equivalently, using the A notation instead of the square notation, one obtains 



(5.4) 







yi,-^{x\y^\---\y'^)+^{x\y^\- 

m 

+ Y,y'^- [2A(x|2/^|---|WT' 

[A {x\y^\ 



'xx\ 



1 2/") + 



m m 



jyhyh 



h = l l2=l 



m m 



-2 5/^A(xy'^|yi| 



L=l ln=l l..=l 



Unfortunately, the equations (5.3) are not solved with respect to the y^^, because in its 
last line, we notice that the y^^^ are mixed with the y^^ . Consequently, we have to solve 
a linear system of m equations with the unknowns yi^ of the form 
(5.5) 



= A^ 



h=i 



«2 = 1 



-X.J-, + X.i. 
y ^ i;'2 w'l y ■ 



for j = 1, . . . , m, where is an abbreviation for the terms appearing in the lines 5, 6, 
7, 8, 9 and 10 of (5.3), or even more compactly, changing the index j to the index k 



(5.6) 



fox k = 1, . . . , m, where is an abbreviation for the terms in the brackets in (5.5). 

Thanks to the assumption that the determinant (3.2) is the identity determinant at 
(x, y) = (0, 0), we deduce that the determinant of the m x m matrix (-Bfi)i<fj<m is 
also the identity determinant at (x, y, y^.) = (0, 0, 0). It follows that the determinant of 
the m X m matrix is nonvanishing in a neighborhood of the origin in the 

first order jet space. Consequently, we can apply the rule of Cramer to solve the y^^ 
explicitely interms of the A'' and of the as follows 



(5.7) 



Bl ■ 


■ A' ■ 


■ Bl 


B"^ ■ 


■ A™ ■ 


■ -Dm 


Bl ■ 


■ B' ■ 


■ Bl 


5f ■ 


■ B"" ■ 


m 



where on the numerator, the only modification of the determinant of the matrix 
(-^fi)i|fi<m is the replacement of its j-th column by the column vector A. We have to 
show that after replacing the A'' and the fifj by their complete expressions, one indeed 
obtains the desired equation (5.4). As in (3.43), we shall introduce a notation for the 
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two m X m determinants appearing in (5.6): we write this quotient under the form 



(5.8) <^ y: 



J 

XX 











Rfc II 
Bm\\ 










II 

m II 



where it is understood that Bl, . . . , . . . , B^ and A are column vectors whose 
index k (for their lines) varies from 1 to m. This representation of determinants em- 
phasizing only its columns will be appropriate for later manipulations. 

Our first task is to compute the determinant in the denominator of (5.8). Recalling 
that we make the notational identification = x, it will be convenient to reexpress 
the Bf in a slightly compacter form, using the total differentiation operator D: 



(5.9) 



■ DX - Xy:, ■ DY' 



Lemma 5.10. We have the following expression for the determinant of the matrix 



(5.11) 



\Yy\ ■ DX - Xyi ■ DY^\ ■ ■ ■ |Fjf;. ■ DX - Xy^ ■ DY''\ 



[DX]""-^ ■ A {x\y^' '-'"^ 



Proof. By multilinearity, we may develope the determinant written in the first line 
of (5.11). Since it contains two terms in each columns, we should obtain a sum of 
2™ determinants. However, since the obtained determinants vanish as soon as the 
column DY'' (multiplied by various factors Xyi) appears at least two different places, 
it remains only (m + 1) nonvanishing determinants, those for which the column DY'' 
appears at most once: 
(5.12) 



\Y^, ■ DX - Xyi ■ DY'^l ■ ■ ■ \Yyi. ■ DX - Xyn. ■ DY^W = 

= [DXr ■ \\Yy\\ ■ ■ ■ \Yyt\\ - [DXr-' Xyi ■ II W^li;' I ■ ■ ■ \Yyt 



- [DXr~^Xy^ ■ \\Yy\\ ■ ■ ■ |yt_i|DF'- 



To establish the desired expression appearing in the second line of (5.11), we factor 
out by [DX]^^^ and we develope all the remaining total differentiation operators D. 
Since = x, we have y° = 1, and this enables us to contract X^ + J2h ^y'l 
Ym^=o yx Xyh . So, we achieve the following computation (further explanations and 
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comments just afterwards): 
(5.13) 



[DX] 



m—l 



Ji=0 

m 

m 



/i=0 



I 



IF* 



/i=0 



1 yl\ - ■ ■ \1 ym-1 I -I 



[DX] 



m—l 



E ^y'l ■ II -^jJ I ■ ■ ■ i^j/" II ~ 



Ji=0 



iJx ' W y^\ y'^ \ ' ' ' \ y"^\\ 

—y^ym ■ II -f yl I ■ ■ ■ I -r ym-1 1 -fa; II ~ Vx ■ II -f I 

= [DXr-'-{X,,-\\Yy\\---\Yyt 

y II I W^W \ 

— y^ym ■ II -T yl I ■ ■ ■ I -r ym-1 II I 

[ =[DXr~'-{Aix\y'\---\yn}. 



Xy..\\Y^'\Y^,\...m 



IK 



iT^fc iT^fc in 

I J- ym-~l I ym II J 



For the passage to the equality of line 4, using the fact that a determinant having two 
identical columns vanishes, we observe that in each of the m sums Yldl=a appearing 
in lines 2 and 3 (including the cdots), there remains only two non- vanishing de- 
terminants. For the passage to the equality of line 7, we just sum up all the linear 
combinations of determinants appearing in lines 4, 5 and 6. Finally, for the passage 
to the equality of line 9, we recognize the development of the fundamental Jacobian 
determinant (3.2) along its first line (X^,, Xyi, . . . , Xym), modulo some permutations 
of columns in the m x m minors. The proof is complete. □ 



Our second task, similar but computationnally more heavy, is to compute the deter- 
minant in the numerator of (5.8). First of all, we have to re-express the A'' defined 
implicitely between (5.3) and (5.5) using the total differentiation operator to contract 
them as follows 
(5.14) 

A' = DX- Y^, - DY' ■X,, + 2J2 Vx " 

* mm 
h = l l2 = l 



DX-Y^i, -DY^ 

xy^ 



X^yll 



DX ■ Yyllyh — DY^ ■ Xyl^yl2 



Replacing this expression of A'' in (5.8), taking account of the expression of 
the denominator already obtained in the second line of (5.11) and abbreviating 
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A(x|y^| • ■ ■ as A, we may write (5.8) in length and then develope it by linear- 
ity as follows 
(5.15) 

\Y^^ ■ DX -Xyi ■ DY''\ ■■■ 
■■■\^ DX- Y^^ - DY'^ ■ Y^^+ 



-1 



-1 



m 



11=1 

m m 



+ 



E y'^ y'^ ■ [^^ ■ K^w^ - ■ ^v^^ 

ll = l l2=l 



yh 



■ ■ ■ \ Yym ■ DX — Xy-m ■ DY 



\Y^^ - DX -Xyl -DY^l--- 



■■■I'DX-Yl-DY'^-X^J--- 



\Yyt-DX-Xym-DY''\\ + 



h=i 

■■■\^ DX ■ Y^ I, - DY^ ■ X„i, I • • • 

■■■\Yyt-DX-Xy,n-DY''\\ + 

m m 

+ E E y'^ ■\K-Dx-Xy.-DY'^\.-. 

■■■\^DX-Y},^,,-DY'-Xy^,y^,\--- 

■ ■ ■ \ Yym ■ DX — Xym • DY^^ 



As it is delicate to read, let us say that lines 2, 3 and 4 just express the j-th colum p A'^\ 
of the determinant \Bik\ ■ ■ A^\ ■ ■ ■ after replacement of A'^ by its complete 

expression (5.14). 

In lines 6, 7, 8; in lines 9, 10, 11; and in lines 12, 13, 14, there are three families 
of m X m determinants containing a linear combination (soustraction) having exactly 
two terms in each column. As in the proof of Lemma 5.10, by multilinarity, we have 
to develope each such determinant. In principle, for each development, we should 
get 2™ terms, but since the obtained determinants vanish as soon as the column DY^ 
(modulo a multiplication by some factor) appears at least twice, it remains only (m+ 1) 
nonvanishing determinants, those for which the column DY^ appears at most once. In 
addition, for each of the obtained determinant, the factor [DXY^~^ appears (sometimes 
even the factor [Z^X]™), so that this factor compensates the factor [DX]™-^ in the 
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numerator. In sum, the continuation of the huge computation yields: 
(5.16) 

' flY . \\V^ I . . . P I . . . IV^ II — 

W y''-' I I II 



-Xyl ■ \\DY''\--- 


1 xa; 


1 IV^ II 
1 ■ ■ ■ Mj/'^ll 








Y . II 1 . . . 

^XX J-yl 11 


DY'' 


1 IV^ 

1 ■ ■ ■ l-fy^^ll 








— Xym ■ \\Yyl \ ■ ■ ■ 


i yk 1 










m. 

ii=i 


Vk 1 


. . |i 1 . . . \yk 

1 zj/'i 1 1 ^ y'" 








m 

h=l 




1 |i 1 Iv'^ 

1 1 -'x-y'il l-^i/' 








Tl 






m 


II Vk 

II V 






— 


m 

-2Y,yxXy-- 

U —1 


VI 


■ ■■p■y^J ■■■IDF 




+ 




m m 


yk 1 . . . |i yk 1 . 

1 y''^y''2 ' 




1 V'fc 




m m 

il = l ^2 = 1 


1 1 y'ly'2l 




lyk 

■ \ Iym 




m m 

"EE 

il = l «2 = 1 


■ y^*!] ■■■p'DF^i 




Wk 

■ \ Iym 




m m 

-EE^^'^^'^^-- 


1 . . . P y'' 1 . 

1 -'^1^2 1 




\DY^ 





To establish the desired expression (5.4), we must develope all the total differentiation 
operators D of the terms DX placed as factor and of the terms DY^ placed in various 
columns of determinants. We notice that in developing DY^, we obtain columns Y'] 

(multiplied by the factor |/^) and for only three (or two) values of / = 0, 1, . . . , m, 
this column does not already appear in the corresponding determinant, so that (m — 
1) determinants vanish and only 3 (or 2) remain nonzero. Taking account of these 
simplifications, we have the continuation 



(5.17) -i/^^^ A = 1 + 11 + 111, 



where the term I is the development of lines 1, 2, 3, 4 of (5.16); the term II is the 
development of lines 5, 6, 7, 8 of (5.16); and the term III is the development of lines 9, 
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10, 11, 12 of (5.16). So we get firstly (further explanations follows): 

m 

^ y Vx Xyl ■ \Yyl I ■ ■ ■ P ■ ■ ■ \Yym || — 



-Xy^-\\yx\--- 

^ , 1 \^ 1 1 ^/^ ^ 1 
- VxXy^ ■ II Vl 


\j 1 . . . 

1 XX \ 

1 J vk 1 

I XX I 


\Yym 

II 

\1 ym II 


~ ViXyJ ■ \\Yyj \ ■ 


■ ■ \^ Y^ 1 

1 xx\ 


|l^m 1 


~ Xxx ■ \\Yyi ■ ■ ■ 


|i y^l . . . 

1 X \ 


\Yym II — 


~ yiXxx ■ \\Yyl\ 




1 V^m II 


— Xym ■ Yyi ■ ■ 


|i 1 . . 

1 xx\ 


■i>;1l- 


Vx ^y"^ Pj/i 


. . . \i 

1 XX 


1 ■ ■ ■ \ Yym 


~ yiXym ■ \\Yyi\ 


. . . u' y'^^ 1 

1 xx\ 





and secondly (we discuss afterwards the annihilation of the underlined terms): 



m m 



(5.19) 



2EE^i^^i^.- yy'\---nxy'^\---K- 

1-,=! 1=0 



ii=i 

m 



h=i 



II 



m 

9 \^ ih Y , . \\V^ I . . . I-? . . . II 

/ ^ ilx ^xyn ||-'j/i| I a; I I J/™ || 

h=l 
m 

9 ih Y , . Ily'^ I . . . \j I . . . ly* 

^ / Ux Ux^xyh II 2/^1 I y^\ \ y 

h=l 

m 

2j2yi^Xy,.. Yy\\...nx,A---\y^ 



h=i 



2j2yxyTXyr.- 


^y" 




\j yfc 

1 xy'i 


1 ■ ■ ■ \ Yym 


7 


m 


Vk 1 




1 







«1=1 



CHARACTERIZATION OF NEWTONIAN FREE PARTICLES 

and where thirdly (we are nearly the end of the proof): 



71 



(5.20) 



m m m 



■yk I I j yk I \yk 



h=l l2=l 1=0 



m m 



yk\ 

X I 



yhyh I 



J2J2y'^'y'^"y'^^y'- 


Vk 1 


\i yk 1 \yk 
1 ^yhyhl l-'j/™ 




ll=l l2 = l 






- S 



Zl=l (2=1 



j Vk I ... I V''^' 

^yhyhl I -'a™ 



m m 



-J2Y.y'^^y'^"^yw--\K 

Zl=l (2=1 



. |i v^l . . . I v^' II — 

I ^ X \ I i/™ II 

E E Xyhyh ■ \\Yyl \ ■ ■ ■ \^ Yyj \ ■ ■ ■ \ Yym 



m m 



h=l l2 = l 
m m 



E E y'^'^y- ■ ^fi I ■ ■ ■ I' ^^1.'^ I ■ ■ ■ 

ll=l l2=l 



m m 



E E y'^' y'^ y^^y- ■ 

Zl=l l2=l 






|i yk 

1 ^yhyl2 


1 ■ ■ ■ \ Yym 




m m 

E E y'i y'^ y^^y- ■ 






j yk 1 







h=l l2 = l 

Now, we explain the annihilation of the underlined terms. Consider I: in the first sum 
Ym^O' terms except only the two corresponding to / = and to / = j are 

annihilated by the other terms with Q] appended: indeed, one must take account of the 
fact that in the expression of I, we have two sums represented by some cdots, the 
nature of which was defined without ambiguity in the passage from (5.15) to (5.16). 
Similar simplifications occur for II and for III. Consequently, we obtain firstly: 



(5.21) I := 



X.-\\Yy\ 


. . . |i y'' I . . . W'^ II -1- 

1 ^ xx\ \ y"^\\ ' 


+ yi Xyj 


Wvk 1 
■ II ^y' 1 


. . . |J yk 1 . . . \yk _ 

\ XX \ \ y™ II 


- Xyl ■ 


Y^\ 

X 1 


|i V^ 1 . . . W'^ 

\ ^ xx\ 1 3/™ II 


- yi Xyj 


II \ 
■ Vy'\ 


. . . |i V'^ 1 . . . 1 V'"" II 
1 ^ xx\ 1 y"Ml 


Xxx ' 1 




\j yk\ \yk 

1 ^ X \ 1 3/™ II 


- yi 


■K\ 


\jyk\ \yk II 
■ ■ ■ 1 1 ■ ■ ■ 1 J/"* II 




\ vk 1 


.\i yk \ . . . \yk\\ _ 

1 ^ xx\ 1 2^ II 


-yiXym 


II vk 

■ II V 


. . . |i yk 1 . . . \yk II . 

1 ^ xx\ \ y^W 1 
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just above, the first two lines consist of the two terms in the sum underlined at the first 
line of (5.18) which are not annihilated; secondly we obtain: 



(5.22) 



+ 



/i=i 



h=i 



II : = 



/i=i 

m 



A^Y^ I . . . 

/i=i 

m 

2 ^ Vx^xyh ■ II^Ji I ■ ■ ■ P^af I ■ ■ ■ l^i/™ || 



li=l 



-2j2yi^yiX,y.,-\\Y^\\---n,'\---\Yi, 



h=i 



m 



h=l 



similarly, the first two lines above consist of the two terms in the sum underlined at the 
first line of (5.19) which are not annihilated; and thirdly we obtain: 



m m 



+ 



+ E E y'^ y^^ y^^ '^y^ ' '^y^^" ' P ^!/'iy'2 1 ■ ■ ■ 1^?^ 



ll=l l2 = l 
m m 



(5.23) III : = 



h=l l2 = l 



■ ■ ■ P' I 

-' X I I -' I 



\Y m 



m m 



- E E y'^ y'^ ^y^ ■ i ■ ■ ■ i' ^i^y^^ i ■ ■ ■ 

/l=l «2 = 1 

m m 

~ ^ ^ ^ ] Xyhyh ■ 11^1 I ■ ■ ■ P ^"1 ■ ■ ■ \ Yyni II - 

Zl=l l2 = l 

m m 

~ E E ^2:^ Xyllyl2 ■ \\Yyl \ ■ ■ ■ \^ Yyj \ ■ ■ ■ \Yym 

ll=l l2 = l 
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I 



\j yk 
I ^yhyh 



h = l l2 = l 

m m 

Collecting the odd lines of (5.21), we obtain exactly (m + 1) terms which correspond to 



the development of the determinant A(x| 



xx\ 



\y^) along its first line, modulo 



permutations of columns of the associated m x m minors; collecting the even lines 
of (5.21), we obtain exactly (m + 1) terms which correspond to the development of 
the determinant —yi ■ A(xx\y^\ ■ ■ ■ |y™) along its first lines, modulo permutations of 
columns of the associated m x m minors. Similar observations hold about II and III. 
In , we may rewrite the final expressions of these three terms: firstly 



(5.24) 



I 

II 



A(x| ■ ■ ■ p xx| ■ ■ ■ ly") - yi ■ A{xx\y^\ ■ ■ ■ ly"), 

m 

2^yi^-A(x|---Px/|---m- 



ii=i 



m m 
'1 = 1 l2 = l 



ni= $^$^yi^yi^-A( 



x\ 



\yn- 



m m 



-yi}2z2y'^'y'^'-^(y''y''\y'\---\y'"^- 

ll=l «2 = 1 

Coming back to (5.17), we obtain the desired expression (5.4). 

The proof of the — technical, though involving only linear algebra — Lemma 3.32 
is complete. □ 



References 

[BK1989] Bluman, G.W.; Kumei, S.: Symmetries and differential equations. Springer- Verlag, 
New- York, 1989. 

[Cal924] Cartan, E.: Sur les varietes a connexion projective. Bull. Soc. Math. France 52 (1924), 
205-241. 

[Chl939] Chern, S.-S.: Sur la geometric d'un systeme d' equations differentielles du second ordre. 
Bull. Sci. Math. 63 (1939), 206-212. 

[CMS1996] Crampin, M.; Martinez, E.; Sarlet, W.: Linear connections for systems of second- 
order ordinary differential equations, Ann. Inst. H. Poincare Phys. Theor 65 (1996), no. 2, 
223-249. 

[Do2000] DOUBROV, B.; Contact invariants of ordinary differential equations, RIMS Kokyuroku 
1150 (2000), 105-113. 

[DNP2005] Dridi, R.; Neut, S.; Petitot, M.: Elie Cartan's geometrical vision or how to avoid 

expression swell, arxiv . org/ abs/math . DG/ 05042 03. 
[Fel995] Pels, M.: The equivalence problem for systems of second-order ordinary differential 

equations, Proc. London Math. Soc. 71 (1995), no. 2, 221-240. 
[G1989] Gardner, R.B.: The method of equivalence and its applications, CBMS-NSF Regional 

Conference Series in Applied Mathematics 58 (SIAM, Philadelphia, 1989), 127 pp. 
[GM2003] Gaussier, H.; Merker, J.: Symmetries of partial differential equations, J. Korean Math. 

Soc. 40 (2003), no. 3,517-561. 



74 



JOEL MERKER 



[GG1983] Gonzalez Gascon, R; Gonzalez Lopez, A.: Symmetries of differential equations, 
IV. J. Math. Phys. 24 (1983), 2006-2021. 

[GL1988] Gonzalez Lopez, A.; Symmetries of linear systems of second order differential equa- 
tions, J. Math. Phys. 29 (1988), 1097-1105. 

[GK01992] Gonzalez Lopez, A.; Kamran, N.; Olver, P.J.: Lie algebras of vector fields in the 
real plane, Proc. London Math. Soc. 64 (1992), 339-368. 

[GTW1989] Grissom, C; Thompson, G.; Wilkens, G.: Linearization of second order ordinary 
differential equations via Cartan's equivalence method, J. Diff. Eq. 77 (1989), no. 1, 1- 
15. 

[Gr2000] Grossman, D.A.: Torsion-free path geometries and integrable second order ODE sys- 
tems, Selecta Math. (N.S.) 6 (2000), no. 4, 399-442. 

[Hal937] Hachtroudi, M.: Les espaces d'elements a connexion projective normale, ActuaUtes 
Scientifiques et Industrielles, 565, Paris, Hermann, 1937. 

[H2001] Hawkins, T.: Emergence of the theory of Lie groups. Springer- Verlag, Berlin, 2001. 

[HK1989] Hsu, L.; Kamran, N.: Classification of second order ordinary differential equations 
admitting Lie groups of fibre-preserving point symmetries, Proc. London Math. Soc. 58 
(1989), no. 3,387^16. 

[IB 1992] Ibragimov, N.H.: Group analysis of ordinary differential equations and the invariance 
principle in mathematical physics, Russian Math. Surveys 47:4 (1992), 89-156. 

[IB1999] Ibragimov, N.H.: Elementary Lie group analysis and ordinary differential equations. 
Mathematical methods in practice, John Wiley & Sons, Chichester, 1999, xviii+347 pp. 

[Lel980] Leach, P.G.L.: Sl(3, R) and the repulsive oscillator, J. Phys. A 13 (1980), 1991-2000. 

[Liel880] Lie, S.: Theorie der Transformationsgruppen, Math. Ann. 16 (1880), 441-528; trans- 
lated in EngUsh and commented in: AcKERMAN, M.; Hermann, R.: Sophus Lie's 1880 
Transformation Group paper. Math. Sci. Press, Brookline, Mass., 1975. 

[Liel883] Lie, S.: Klassifikation und Integration vo gewohnlichen Differentialgleichungen zwischen 
X, y, die eine Gruppe von Transformationen gestaten I-IV. In: Gesammelte Abhandlungen, 
Vol. 5, B.G. Teubner, Leipzig, 1924, pp. 240-310; 362-427, 432^48. 

[MS2001] Mahomed, F.M.; Son, C.W.: Linearization criteria for a system of second-order differ- 
ential equations, Internat. J. Non-Linear Mech. 36 (2001), no. 4, 611-611. 

[Ma2003] Mardare, S.: On isometric immersions of a Riemannian space under a weak regularity 
assumption, C. R. Acad. Sci. Paris, Ser. I 337 (2003), 785-790. 

[M2004] Merker, J.: Explicit differential characterization of PDE systems pointwise equivalent 
to Yxjixi2 = 0, 1 ^ ji,j2 ^ n ^ 2, arxiv.org/abs/math.DG/0411637. 

[N2003] Neut, S.: Implantation et nouvelles applications de la methode d' equivalence d'Elie 
Cartan, These, Universite Lille 1, October 2003. 

[NS2003] NUROWSKY, P.; SPARLING, G.A.J.: 3-dimensional Cauchy-Riemann structures and 2'^'^ 
order ordinary differential equations, e-print arXiv : math .DG/0306332j 

[011986] Olver, P.J.: Applications of Lie groups to differential equations. Springer Verlag, New 
York, 1986. xxviH-497pp. 

[OL1995] Olver, P.J.: Equivalence, Invariance and Symmetries. Cambridge University Press, Cam- 
bridge, 1995, xviH-525 pp. 
[Stel982] Sternberg, S.: Differential geometry. Chelsea, New York, 1982. 

[Stk2000] Storm ARK, O.: Lie 's structural approach to PDE systems. Encyclopaedia of mathematics 
and its applications, vol. 80, Cambridge University Press, Cambridge, 2000, xvH-572 pp. 

[TR1896] TreSSE, A.: Determination des invariants ponctuels de I'equation differentielle du second 
ordre y" = y, y'), Hirzel, Leipzig, 1896. 



CNRS, Universite de Provence, LATP, UMR 6632, CMI, 39 rue Joliot-Curie, F- 13453 
Marseille Cedex 13, France. , Inte met: \http : / /www . cmi . univ-mrs . fr /^merker /index . html| 
E-mail address: merker@cmi . univ-mrs . f r 



